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FOREWORD 


The  Symposium  on  Aeroelastic  and  Dynamic  Modeling  Technology,  sponsored  by  the 
Air  Force  Flight  Dynamics  Laboratory  of  the  Research  and  Technology  Division.  Air 
Force  Systems  Command,  in  association  with  the  Dynamics  and  Aeroelasticity  Research 
Panel  of  the  Aerospace  Industries  Association,  was  held  in  Dayton,  Ohio,  on  23-24-25 
September  1963.  The  purposes  of  the  Symposium  were  to:  (a)  review  the  state-of-the-art, 

(b)  present  the  latest  developments  and  applications  for  solution  of  current  problems, 

(c)  determine  the  maximum  extent  to  which  modeling  techniques  can  be  employed  for 
efficient  design  and  substantiation  of  full  scale  vehicle  safety  and  reliability  and  (d)  pro¬ 
vide  a  basis  for  future  applied  research  for  aeroelastic  and  dynamic  modeling  technology. 
The  symposium  subject  matter  covers  the  theory,  design  and  testing  of  both  static  and 
dynamic  elastic  models  used  in  the  development  of  aircraft,  missiles,  aerospace  vehicles 
and  equipment. 

The  proceedings  are  arranged  in  two  parts:  Part  I  contains  the  unclassified  papers 
grouped  by  subject  as  presented,  and  Part  II  contains  the  classified  papers  and  is  Confi¬ 
dential.  The  authors’  names  and  company  affiliation  are  given  at  the  beginning  of  each 
paper. 

An  effort  has  been  made  to  minimize  errors  and  omissions,  however,  if  any  serious 
discrepancies  are  noted,  they  should  be  Drought  to  the  attention  of  Mr.  M.  H.  Shirk,  FDDS, 
AF  Flight  Dynamics  Laboratory,  WPAFB,  Ohio.  Necessary  errata  sheets  will  be  prepared 
and  distributed. 

The  success  of  the  Symposium  in  attaining  its  purposes  and  the  subsequent  compilation 
of  this  report  for  wide  distribution  as  a  reference  document  were  made  possible  by  the  fine 
contribution  of  the  Symposium  speakers,  session  chairmen  and  authors  of  the  technical 
papers.  The  AFFDL  and  the  AIA  gratefully  acknowledge  these  contributions  and  express 
their  appreciation. 

This  report  is  classified  CONFIDENTIAL  because  it  contains  information  from  which 
performance  of  military  weapon  systems  under  development  can  be  obtained,  the  unauthor¬ 
ized  disclosure  of  which  is  considered  prejudicial  to  the  defense  interests  of  the  United 
States. 
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ABSTRACT 


■Hie  Confidential  proceedings  of  the  Air  Force  Flight  Dynamics 
Laboratory  and  Aerospace  Industries  Association  Symposium  on  Aero- 
elastic  and  dynamic  Modeling  Technology  are  presented  in  this  report. 

Modem  day  aerospacecraft  requirements  are  calling  more  and 
more  for  systems  which  are  greatly  influenced  by  aero  and  thermo  - 
elastic  effects.  In  line  with  this  trend,  model  testing  technology 
has  been  required  to  keep  pace  by  advancing  the  state-of-the-art 
in  the  simulation  of  elastic,  thermal,  aerodynamic  and  dynamic 
effects.  The  objective  of  this  symposium  was  to  bring  together  the 
members  of  the  scientific  community  who  are  intimately  concerned  with 
this  area  of  endeavor,  thereby  promoting  a  mutual  understanding  of 
current  developments  and  problem  areas.  The  papers  presented  at  the 
symposium  were  divided  into  appropriate  technical  area  sub-groups 
and  individual  sessions  were  devoted  to  each.  These  sub-groups  were: 
Theory  and  Design,  Model  Testing  Techniques,  Dynamic  Loads  and  Aero- 
elastic  Applications,  and  Structural  Design  Applications  with  two 
classified  sessions  on  Aerospace  Vehicle  Applications  and  Aircraft 
Applications . 


PUBLICATION  REVIEW 


The  publication  of  this  report  does  not  constitute  approval  by 
the  Air  Force  of  the  findings  or  conclusions  contained  herein.  It 
is  published  only  for  the  exchange  and  stimulation  of  ideas. 


WILLIAM  C.  NIELSEy 
Colonel,  USAF 
Symposium  Chairman 
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DYNAMIC  STABILITY  EQUIPMENT  AND  TECHNIQUES 
FOR  TESTING  AT  HYPERSONIC  SPEEDS* 

By 

L.  K.  Ward  and  R.  H.  Urban 
von  K^rm^n  Gas  Dynamics  Facility 
ARO,  Inc. 

ABSTRACT 


Dynamic  stability  balances,  both  free  and  forced  oscil¬ 
lation,  developed  at  the  VKF  specifically  for  hypersonic 
measurements  are  described  and  discussed.  Problems  encoun¬ 
tered  in  measuring  hypersonic  model  damping,  and  their  solu¬ 
tions,  are  mentioned  in  addition  to  defining  areas  in  data 
reduction  procedures  where  caution  should  be  exercised.  Data 
obtained  on  several  re-entry  configurations,  primarily  at 
Mach  number  10,  are  reviewed,  pointing  out  trends  found  with 
variation  of  Reynolds  number,  model  frequency  of  oscilla¬ 
tion,  model  amplitude  of  oscillation,  and  model  geometry. 
Finally,  areas  of  interest  for  future  dynamic  stability  re¬ 
search  programs  are  discussed. 
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NOMENCLATURE 

A 

Reference  area  (ird2/4),  ft2 

Cm 

Pitching  moment/q^Ad 

cmq 

dCm/d (qd/2VOT) ,  damping-in-pitch  derivative,  1/rad 

Cma 

Slope  of  the  pitching-moment  curve,  1/rad 

C“d 

dCm/d(ad/2Voo) ,  damping-in-pitch  derivative,  1/rad 

i 

C7R 

Cycles  to  damp  to  amplitude  ratio  R,  cycles 

* 

d 

Reference  length  (noted  model  diameter),  ft 
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f  Frequency,  cycles/sec 

I  Mass  moment  of  inertia,  slug-ft2 

in  Natural  logarithm 

M  Free-stream  Mach  number 

00 

M0  Angular  restoring-moment  parameter,  ft-lb/rad 

/ 

M^ ,  M^  Angular  viscous-damping  parameters,  f t-lb-sec/rad 
q  Angular  velocity 

q  Free-stream  dynamic  pressure,  lb/ft2 

00 

R  Ratio  of  the  amplitude  of  a  damped  oscillation  at  a 

given  time  to  the  initial  amplitude 

Re  Reynolds  number,  based  on  model  reference  length,  d 

rb  Model  base  radius,  in. 

rQ  Model  nose  radius,  in. 

T  Input  torque,  ft-lb 

t  Time,  sec 

V  Free-stream  velocity,  ft/sec 

00 

a  Angle  of  attack,  rad 

a  Time  rate  of  change  of  angle  of  attack,  rad/sec 

0  Angular  displacement,  rad  or  deg 

8  Angular  velocity,  rad/sec 

0  Angular  acceleration,  rad/sec2 

cjd  Angular  frequency,  rad/sec 

SUBSCRIPTS 


a 

Atmospheric  conditions 

o 

Maximum  conditions 

V 

Vacuum  conditions 

w 

Wind-on  conditions 
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1.0  INTRODUCTION 


The  need  to  simulate  the  dynamic  motion  of  proposed 
flight  and  re-entry  vehicles  in  the  supersonic  and  hypersonic 
speed  regimes  and  thus  provide  data  to  properly  evaluate  the 
damping  characteristics  of  such  vehicles  has  led  to  the  de¬ 
velopment  of  several  dynamic  stability  balances  at  the  von 
K^rm^n  Gas  Dynamics  Facility  (VKF)  ,  Arnold  Engineering  De¬ 
velopment  Center  (AEDC),  Air  Force  Systems  Command  (AFSC)  for 
use  in  the  supersonic  and  hypersonic  wind  tunnel  facilities. 

Initial  development  provided  sting-supported,  low  ampli¬ 
tude,  forced  oscillation  balances  for  use  in  the  supersonic 
wind  tunnels.  Modification  of  this  equipment  by  the  addition 
of  air  and  water  cooling  allowed  testing  in  the  hypersonic 
tunnels.  Parallelling  this  hypersonic  capability  was  the 
development  of  several  free  oscillation  balances  which  pro¬ 
vided  model  oscillation  amplitudes  ranging  up  to  ±90  deg. 
Problems  encountered  with  balance  tare  damping  in  the  hyper¬ 
sonic  speed  regime  prompted  development  of  a  high  amplitude 
forced  oscillation  balance  and  a  gas  bearing  pivot  for  free 
oscillation  measurements. 

Dynamic  stability  data  have  been  successfully  obtained 
for  numerous  test  configurations  in  the  Mach  number  range 
from  1.5  to  10  in  tests  sponsored  by  the  Ballistic  Systems 
Division  (BSD)  and  the  Aeronautical  Systems  Division  (ASD)  of 
the  AFSC,  the  Army  Ballistic  Missile  Agency  (ABMA) ,  now  known 
as  the  Army  Missile  Command,  and  the  National  Aeronautics  and 
Space  Administrations  Manned  Spacecraft  Center  (NASA-MSC). 
Results  of  some  of  these  test  programs  will  be  presented. 

Dynamic  damping  data  obtained  in  the  hypersonic  speed 
regime  have  indicated  several  areas  of  disagreement  with  the 
predicitions  of  the  Newtonian  impact  theory  (Ref.  1).  These 
data  are  presented  and  briefly  discussed  along  with  work  pro¬ 
grams  now  in  progress  at  the  VKF. 


2.0  APPARATUS 


2.1  WIND  TUNNELS 

The  40-in.  supersonic  and  50-in.  hypersonic  wind  tunnels 
are  closed-circuit,  variable  density  wind  tunnels,  whereas 
the  12-in.  supersonic  and  12-in,  hypersonic  wind  tunnels  are 
of  the  intermittent,  variable  density  type. 
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The  40-Inch  Supersonic  Tunnel  (A)  (Fig.  1)  has  a  flexi- 
ble-plate-type  nozzle  which  is  automatically  controlled  to 
produce  Mach  numbers  from  1.5  to  6.  The  tunnel  operates  at 
maximum  stagnation  pressures  ranging  from  29  to  200  psia  at 
Ma,  *■  1-5  to  6,  respectively,  and  at  stagnation  temperatures 
up  to  300°F  (Moo  “  6)  .  Minimum  operating  pressures  are  about 
one- tenth  of  the  maximum. 

The  50-Inch  Mach  8  Tunnel  (B)  (Fig.  2)  and  the  50-Inch 
Mach  10  Tunnel  (C)  (Fig.  3)  have  contoured  axisymmetric  noz¬ 
zles.  Tunnel  B  operates  at  stagnation  pressures  ranging 
from  100  to  900  psia  and  at  stagnation  temperatures  up  to 
900°F.  Tunnel  C  operates  at  stagnation  pressures  ranging 
from  200  to  2000  psia  and  at  stagnation  temperatures  up  to 
1450°F . 

The  12-Inch  Supersonic  Tunnel  (D)  (Fig.  4)  has  a  manu¬ 
ally  adjusted,  flexible-plate-type  nozzle.  The  tunnel  oper¬ 
ates  at  Mach  numbers  from  1.5  to  5,  at  stagnation  pressures 
from  about  5  to  60  psia,  and  at  stagnation  temperatures  of 
about  80°F . 

The  12-Inch  Hypersonic  Tunnel  (E)  (Fig.  5)  has  a  nozzle 
formed  by  manually  adjusted  contoured  throat  blocks  and 
flexible  plates.  The  tunnel  operates  at  Mach  numbers  from 
5  to  8  at  maximum  stagnation  pressures  from  400  to  1600  psia, 
respectively,  and  at  stagnation  temperatures  up  to  about 
940°F. 

A  complete  description  of  the  tunnels  and  airflow  cali¬ 
bration  information  may  be  found  in  Ref.  2. 


2.2  DYNAMIC  BALANCES 

Several  one-degree-of-freedom  dynamic  balances  are 
available  at  the  VKF  and  are  summarized  in  Table.  1.  Both 
the  free  oscillation  and  forced  oscillation  techniques  are 
used  for  measuring  the  damping-in-pitch  or  damping- in-yaw 
derivatives . 

2.2.1  Forced  Oscillation 

Figure  6  shows  the  sting-balance  assembly  of  the  low 
amplitude,  forced  oscillation  balance.  In  the  case  of  the 
high  amplitude,  forced  oscillation  balance  the  identical 
sting  support  system  is  used.  Each  system  is  forced  to  os¬ 
cillate  through  linkage  by  an  electromagnetic  shaker  motor 
located  in  the  enlarged  portion  of  the  sting.  The  low  am¬ 
plitude  (±3  deg)  system  (Fig.  7)  uses  a  cross-flexure  pivot 
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which  provides  exceptionally  low  values  of  structural  damp¬ 
ing.  The  angular  displacement  of  the  model  is  measured  by 
the  use  of  a  strain-gage  bridge  mounted  on  a  cross-flexure. 

The  input  torque  to  the  system  is  measured  with  a  strain- 
gage  bridge  mounted  on  the  torque  beam  (Fig.  7);  hence,  the 
damping  torque  contributed  by  the  shaker  motor,  push  rod, 
and  linkage  assembly  is  not  measured. 

The  high  amplitude  (±12  deg) ,  forced  oscillation  bal¬ 
ance  (Fig.  8)  has  an  oscillating  shaft,  supported  by  to 
angular  contact  bearings,  for  the  model  pivot  axis.  The 
bearings  are  spring  loaded  axially  to  eliminate  bearing 
free-play.  The  shaft  is  forced  to  oscillate  by  a  flexure- 
connected  push  rod  which  is  driven  by  the  shaker  motor  pre¬ 
viously  mentioned.  The  push-rod-to-shaft  connection  is 
made  by  means  >f  a  small  flexure,  thus  eliminating  the  use 
of  pin  joints.  Removable  I-beam  section  torque  members  are 
connected  to  the  shaft  on  the  model  side  of  the  bearings  and 
support  the  model  mounting  bulkhead.  Balance  restoring  moments 
are  produced  by  the  two  tapered  beams  which  are  fixed  to  the 
sting  and  connected  to  the  model  mounting  bulkhead  by  two 
sets  of  flexures.  Points  B  are  the  attachment  points  of  one 
set  of  flexures  to  the  bottom  restoring-moment  beam,  and 
points  B'  are  the  attachment  points  of  these  flexures  to  the 
model  mounting  bulkhead.  This  connecting  flexure  arrange¬ 
ment  was  designed  to  minimize  the  high  stress  levels  that 
are  present  in  flexures  of  a  shorter  length  when  forced  to 
bend  through  the  large  displacement  angles  required  in  this 
system.  View  AA  shows  the  balance  displaced  12  deg.  Model 
angular  displacement  is  measured  with  a  strain-gage  bridge 
mounted  on  the  restoring-moment  beam,  and  the  input  torque 
to  the  system  is  measured  by  a  strain-gage  bridge  mounted  on 
the  "I"  section  of  the  torque  member. 

Both  forced  oscillation  balance  systems  are  equipped 
with  an  electronic  feedback  control  system  (Ref.  3)  to  pro¬ 
vide  positive  amplitude  control  for  either  dynamically  sta¬ 
ble  or  unstable  models. 


2.2.2  Free  Oscillation 

Several  free  oscillation  balances  have  been  developed 
at  the  VKF  capable  of  both  low  and  high  amplitude  measure¬ 
ments.  Both  cross-flexures  and  bearings  have  been  used  as 
the  pivot.  Ball  bearings  may  be  used  at  the  lower  Mach  num¬ 
bers  for  configurations  having  a  sizable  amount  of  damping, 
but  their  adverse  damping  characteristics  make  ball  bearings 
unsuitable  for  hypersonic  measurements. 

Because  of  the  need  for  a  more  feasible  free  oscilla¬ 
tion  pivot,  a  gas  bearing  (Ref.  4)  was  developed  at  the  VKF. 
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Figure  9  shows  a  sketch  of  the  bearing,  and  as  may  be  noted, 
the  plenum  is  located  in  the  fixed  portion  of  the  assembly 
and  the  gas  escapes  through  eight  equally  spaced  orifices 
to  "float"  the  outer  ring  to  which  the  model  may  be  attached. 
The  bearing  will  easily  support  a  300-lb  axial  load  at  a 
plenum  pressure  of  300  psig  and  may  be  sting  mounted  or 
mounted  by  a  transverse  rod  support  to  allow  free  oscilla¬ 
tion  measurement  of  model  damping.  The  tare  damping  of  the 
bearing  is  extremely  low,  as  may  be  seen  by  inspection  of 
Fig.  10  which  shows  that  for  a  typical  re-entry  configura¬ 
tion  at  =  10  the  gas  bearing  damping  is  two  orders  of 
magnitude  lower  than  the  tare  damping  of  a  cross-flexure 
pivot.  The  gas  bearing  balance,  when  mounted  on  a  sting 
support,  is  restricted  in  oscillation  amplitude  only  by  the 
model's  physical  size  with  respect  to  the  bearing.  By  sup¬ 
porting  the  bearing  on  a  transverse  rod,  there  are  no  ampli¬ 
tude  restrictions  on  the  model,  as  may  be  seen  in  Fig.  11. 

The  angular  transducer  (Fig.  12)  designed  for  the  bear¬ 
ing  is  similar  to  the  one  used  in  previous  VKF  free-oscilla- 
tion  balances.  The  transducer  simply  consists  of  an  eccen¬ 
tric  mounted  on  the  movable  portion  (outer  ring)  of  the 
bearing  and  the  two  "E"  cores  mounted  180  deg  apart  and  over 
the  eccentric.  As  the  bearing  is  rotated,  the  air  gap  be¬ 
tween  the  eccentric  and  the  coils  in  the  "E"  cores  changes 
and  produces  an  analog  signal  proportional  to  the  angular 
displacement.  No  torque  is  produced  by  the  transducer. 


3 . 0  PROCEDURE 


3.1  FORCED  OSCILLATION 

The  motion  of  a  single-degree-of-freedom,  forced  oscil¬ 
lation  system  may  be  defined  by  the  equation 

10  -  -  M gd  =  T  cos  art  (1) 

It  may  readily  be  shown  that,  for  constant  amplitude 
motion  at  the  undamped  natural  frequency  of  the  model-balance 
system,  the  required  forcing  torque  (T  cos  cut)  is  precisely 
equal  to  the  damping  torque  (M §9) .  The  method  for  computing 
the  dimensionless  damping-in-pitch  derivatives  for  constant 
amplitude,  forced  oscillation  tests  is  as  follows: 

Ma9  =  T  cos  cot  (at  the  undamped  natural 
frequency  of  the  system) 

Ma@  <x>  =  T 
0  o 


(2) 

(3) 
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Mg  =  T/eoo> 

(4) 

Mg'  =  Mgw  -  M ev  •  ay^w 

(5) 

+  Cmd  *  M-  •  Wjqjd* 

(6) 

The  expression  ,for  obtaining  the  aerodynamic  viscous 
damping  parameter  M^  is  based  on  the  fact  that  the  structur¬ 
al  damping  of  a  flexure  pivot  varies  inversely  with  the  fre¬ 
quency  of  oscillation,  as  was  shown  by  Welsh  and  Ward  in 
Ref.  5.  Further  experiments  at  the  VKF  have  revealed  that 
all  pivot  systems  investigated,  having  structural  stiffness, 
have  this  inverse  relationship  with  frequency. 

All  forced  oscillation  balance  systems  are  equipped 
with  transducers  which  provide  continuous  time  resolved  out¬ 
puts  which  are  proportional  to  the  model  oscillation  ampli¬ 
tude  and  input  torque  of  the  system .  The  signals  are  digi¬ 
tized  by  a  Beckman  210  high-speed  digitizer  and  relayed  to 
an  IBM  7070  computer  for  data  reduction. 

Operation  of  balances  in  the  hypersonic  tunnels  requires 
that  balance  and  balance  support  assembly  be  adequately 
cooled  to  withstand  the  necessarily  high  stagnation  tempera¬ 
tures.  Water  jackets  provide  sufficient  cooling  for  the  bal¬ 
ance  supports ;  however ,  water  may  not  be  used  to  cool  the 
balance  since  the  structural  damping  of  the  system  would  be 
affected.  Cooling  of  the  balance  has  been  done  with  the  use 
of  high  pressure  air  sprays  directed  on  the  model  attachment 
bulkhead,  as  indicated  in  Fig.  13.  The  cooling  air  is  turned 
off  during  time  periods  when  data  are  being  obtained.  In  the 
case  of  tunnel  C  the  model  is  retracted  from  the  airflow 
during  the  air  cooling  periods. 

As  indicated  by  Eq.  (5),  structural  damping  values  are 
evaluated  at  vacuum  conditions.  This  is  done  to  eliminate 
the  still-air  contribution  encountered  at  atmospheric  condi¬ 
tions.  Normally,  the  still-air  contributions  may  be  ignored 
when  testing  at  the  lower  Mach  numbers  (M®  <  4)  but  must  be 
considered  when  testing  at  the  higher  Mach  numbers  because 
the  results  will  be  in  error,  as  is  indicated  in  Fig.  14. 


3.2  FREE  OSCILLATION 

The  equation  of  motion  of  a  model  suspended  in  an  air- 
stream  and  restricted  to  one  degree  of  freedom  and  small  am¬ 
plitudes  can  be  written  as 


19  -  Mq0  -  M0e  =  0 


(7) 
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This  equation,  of  course,  is  true  only  for  a  completely 
linear  system  with  viscous  damping.  The  damping  parameter 
for  motion  defined  by  Eq.  ( 7 )  can  be  written  as 

M-  =  (8) 

6  C7R 

From  this  expression  it  is  seen  that,  knowing  the  model  mo¬ 
ment  of  inertia,  it  is  only  necessary  to  determine  the  num¬ 
ber  of  cycles  to  damp  to  a  given  ratio  (Cyp)  and  the  fre¬ 
quency  of  oscillation  (f)  to  evaluate  the' damping  parameter 


In  a  true,  viscous,  damped  system  the  amplitude  of  the 
motion  will  decrease  exponentially,  e.g.,  a  plot  of  the 
natural  logarithm  of  the  envelope  curve  against  cycles  of 
oscillation  will  be  linear.  In  reducing  data  for  models 
oscillating  over  large  amplitudes  and  subject  to  non-line¬ 
arities,  it  is  necessary  to  ensure  that  the  amplitude  range 
can  be  approximated  by  the  linear  expression  in  Eq.  (8). 

In  the  case  of  ball  bearing  damping,  the  procedures  for 
evaluating  their  damping  characteristics  is  given  in  Ref.  6. 
After  the  damping  has  been  evaluated  at  the  proper  frequency 
and  amplitude,  the  subtraction  may  be  made  from  the  wind-on 
value  to  yield  the  aerodynamic  value. 

To  this  poi,nt  it  has  been  assumed  that  the  model  bal¬ 
ance  system  has  no  structural  stiffness.  In  the  case  where 
structural  stiffness  is  present  the  aerodynamic  value  must 
be  determined  as  follows 

M-  =  2UnR[(f/CyR)w  -  (f/CyR)v  •  f/fj  (9) 

This  expression  provides  the  necessary  correction  for  the 
variation  of  structural  damping  with  frequency  as  previously 
discussed . 

Most  free  oscillation  measurements  to  be  made  in  the 
future  will  be  made  using  the  gas  bearing  pivots;  therefore, 
the  balance  tare  damping  will  be  extremely  low  and  in  most 
cases  may  be  ignored. 

As  with  the  forced  oscillation  systems  the  full  data 
reduction  may  be  done  on  the  IBM  7070.  Cooling  of  the  free 
oscillation  equipment  is  also  handled  in  the  manner  previ¬ 
ously  described. 
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4.0  EXPERIMENTAL  RESULTS 


Numerous  dynamic  stability  test  programs  have  been. con¬ 
ducted  at  the  VKF  over  the  past  several  years  on  ballistic 
re-entry  vehicles,  and,  in  general,  results  obtained  at  the 
higher  Mach  numbers  have  agreed  reasonably  well  with  the 
Newtonian  impact  theory  (Ref.  1).  In  addition,  no  appreci¬ 
able  effect  was  found  on  the  damping  derivatives  which  was 
attributable  to  changing  the  reduced  frequency  parameter 
(cDd/^Vaj)  .  Tests  conducted  on  a  model  of  the  Pershing  re¬ 
entry  vehicle  configuration  (Ref.  7)  in  the  12-Inch  Hyper¬ 
sonic  Tunnel  (E)  showed  excellent  agreement  with  impact 
theory,  as  may  be  noted  in  Fig.  15. 

Data  obtained  on  models  of  a  typical,  second  generation 
re-entry  vehicle  in  the  VKF  tunnels  A,  B,  D,  and  E  at  Mach 
numbers  4  through  8  are  presented  in  Fig.  16.  The  reduced 
frequency  parameter,  cDd/2V00,  was  varied  by  changing  both 
model  scale  and  frequency  and  ranged  from  0.0039  to  0.007. 

No  apparent  effect  was  noted.  A  point  of  interest  here  is 
that,  in  order  to  produce  a  dynamically  stable  conf iguration 
at  the  lower  Mach  numbers ,  nose  roughness  was  needed  to  pro¬ 
vide  a  turbulent  boundary  layer  and  eliminate  the  flow  sepa¬ 
ration  on  the  cylinder  caused  by  the  10-deg  flare. 

Investigations  of  three  model  configurations  of  the 
Skybolt  re-entry  vehicle  at  Mach  10  in  the  50-Inch  Mach  10 
Tunnel  (C)  (Ref.  8)  indicated  a  pronounced  Reynolds  number 
effect  (Fig.  17).  During  the  tests  the  distance  from  the 
model  base  forward  to  the  pivot  axis  was  held  constant  while 
the  model  nose  radius  was  varied.  The  impact  theory  for 
Models  A,  B,  and  C  shown  in  the  band  did  predict  a  reduction 
in  stability  with  increasing  nose  radius  although  not  of  the 
magnitude  found  by  experiment.  The  theory  should  become  in¬ 
creasingly  valid  as  Reynolds  number  increases,  yet  the  dis¬ 
agreement  with  the  experimental  data  is  the  greatest  at  the 
higher  Reynolds  numbers. 

In  a  more  recent  test  program,  damping  derivatives  were 
measured  on  a  proposed  re-entry  vehicle  and  variations  caused 
by  Reynolds  number,  angle  of  attack,  and  model  oscillation 
frequency  were  investigated  for  low  oscillation  amplitudes 
(±2  deg)  at  Mach  number  10  in  Tunnel  C.  Data  were  obtained 
at  model  frequencies  ranging  from  about  4  to  18  cycles  per 
second  and  over  the  full  Reynolds  number  range  of  the  tun¬ 
nel.  These  data  are  reported  in  full  in  Ref.  9,  and  the 
zero  angle-of-attack  values  are  presented  in  Fig.  18.  As 
may  be  noted,  the  higher  value  of  Cmq  +  occurs  at  the 
lowest  Reynolds  number  and  the  highest  oscillation  frequency 
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and  is  an  order  of  magnitude  higher  than  data  obtained  at  the 
reverse  of  these  conditions.  These  data  show  that,  for  values 
of  0X1/2V00  £  0.0019,  Reynolds  number  has  no  appreciable  effect 
on  the  damping  derivatives.  As  a>d/2Voo  is  idcreased  above 
0.0019  the  damping  derivatives  increase  and,  in  addition,  in¬ 
crease  with  decreasing  Reynolds  number.  The  dotted  line, 
representing  a  "plane"  corresponding  to  the  impact  theory, 
shows  fair  agreement  at  all  Reynolds  numbers  corresponding 
to  values  of  O)d/2V00  £  0.0019. 

It  should  again  be  pointed  out  that  the  above  variations 
with  the  reduced  frequency  parameter  ( a^d/2Vcx> )  were  made  by 
holding  the  model  scale  (d)  constant  and  the  free-stream  ve¬ 
locity  (VM)  essentially  constant.  Additional  experiments 
should  be  conducted  to  determine  if  the  abovementioned  trends 
with  0X1/2VO3  are  valid  for  conditions  where  model  frequency  is 
held  constant  and  geometry  is  varied. 


5.0  FUTURE  TEST  PROGRAMS 


In  a  continuing  effort  to  provide  necessary  test,  equipment 
for  use  in  obtaining  dynamic  damping  data  at  the  VKF,  several 
internal  development  programs  have  been  initiated.  One  such 
program  proposes  the  measurement  of  model  damping  by  obtain¬ 
ing  continuous,  time-resolved,  model  surface  pressures  on  a 
model  oscillating  at  constant  oscillation  amplitudes.  Data 
will  be  obtained  at  amplitudes  ranging  from  ±2  to  ±10  deg 
with  initial  testing  being  done  in  the  tunnel  A  facility  at 
Mach  4.  Variable  reluctance,  dual  coil,  pressure  transducers, 
as  described  in  Ref.  2,  will  be  used  and  will  be  located  at 
the  model  surface  to  minimize  lag  times. 

The  effects  of  ablation  on  the  aerodynamic  characteris¬ 
tics  of  re-entry  vehicles  has  prompted  the  study  of  a  method 
to  provide  damping  data  using  mass  injection  at  the  model 
surface.  Investigations  are  currently  being  conducted  to 
determine  a  means  of  transporting  the  mass  across  the  balance 
pivot  axis  without  appreciably  effecting  the  tare  damping. 


6 . 0  SUMMARY 


A  relatively  large  selection  of  test  equipment  is  avail¬ 
able  at  the  VKF  for  obtaining  model  damping  data  at  Mach  num¬ 
bers  ranging  from  1.5  through  10.  Both  the  free  oscillation 
and  the  forced  oscillation  techniques  are  employed.  The  mod¬ 
els  may  be  supported  by  either  a  sting  or  a  transverse  rod. 
Experimental  dynamic  stability  data  obtained  at  hypersonic 
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speeds  point  out  the  inadequacies  in  the  present-day  theory 
and  indicate  the  need  for  continued  wind  tunnel  investiga¬ 
tions  in  support  of  the  development  of  proposed  flight  and 
re-entry  vehicles. 
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Summary  of  Balances  Currently  Available  for  Dynamic  Stability  Testing 


CONFIDENTIAL 


c 

o 

§  s* 

l fS 

o 

c S  • 

o 

r-^ 

o 

c  *. 

'w  C 

» 

r— 

O  CSJ 
ro  cm 

■+* 

cz  n 

£  - 
00 

0) 

of 

.. 

irT 

n- 

cn 

ITS 

^r*  csT 

# 

*cr 

CM 

CSI 

Fund 

Model 

Q) 

C  L  n 

g  a>  TO  cT 
■g  c  c  ,= 
S  -  TO  o 

a> 

<V 

c 

</> 

<o 

V_ 

.11 

a> 

> 

a> 

> 

C7 

cn 

C7> 

cn 

cn 

Q-  3 
GO 

C 

_C 

_c 

C 

jr  -a 
2  o 

c 

i  E 

00 

go 

go 

OO 

h-  Qi 

GO 

f—  or 

a> 

cn 

or 

cn 

cn 

cn 

<✓> 

</) 

o 

0. 

G. 

c 

c 

C 

C 

C 

3 

tsi  — T 

■— 

>s 

X 

©  K 

—  TO 

—  TO 

—  TO 

u. 

w- 

TO  <V 

TO  <L> 

to  a> 

TO  <1> 

to  a> 

o 

u. 

O  LL. 

CO  JO 

CO  CO 

co  CO 

O  CO 

O  00 

f  f 

CNJ 

o 

3 

CSJ 

§ 

CO 

CO 

CNJ 

—  clT 

-H 

-H 

-H 

+i 

-H 

-H 

-H 

O. 

P  c 

t 

t 

t 

t 

t 

t 

t 

jz  to 

r— 4 

CO 

LO 

l— H 

LO 

**■  o: 

41 

-H 

+1 

■H 

-H 

-H 

M 

o 

JU 

O 

o 

c 

c 

c 

c 

o 

o 

.  o 

o  o 

o  O 

to 

</> 

JO  .— < 

OO 

LO  f— 1 

\JY  f—i 

&  " 

k_ 

a> 

t- 

Q> 

.C  -C 

-C 

a?  ^ 

a>  -c 

>  = 

CL 

3 

CL 

3 

£  £ 

£ 

§■» 

§-£ 

Q  C 

on 

sz 

.c  .c  .c 

S 

x: 

00  ^ 
JZ  iz 

OO  S 
x:  jc 

o 

o 

O  U 

o 

O  G? 

o  u 

c 

c 

c  c 

c 

c  c 

c  c 

6 

i  i  o 

i  cb 

Ai 

m  lt\ 

LO 

^  LO 

'T  lO 

c 

c 

c 

C 

c 

c 

C 

o> 

o 

_o 

o 

_o 

o 

y  o> 

CL 

*o 

a> 

s  s 

■si 

T3  J2 
(L>  - — 

_TO 

TO 

TO 

TO 

TO  1 

o 

Ll_  u 

S  O 

if  -0 

CO  <-> 

aj  — 

a>  u 

Qi  O 

<u  — 

a>  o? 

L_  ^ 

o 

o  JC 

W  ^ 

L.  JO 

L.  «/» 

LL. 

o  O 

u_  O 

LL.  O 

u-  O 

u_  O 

u.  O 

U_  O 

c 

a>  o 
u 

c  £ 

TO  C 

“  cn 

TO  •— 

a> 

■o 

3 

"a. 

E 

< 

0> 

*o 

3 

"5. 

E 

< 

a> 

■o 

3 

"5. 

E 

C 

a> 

■o 

3 

"5. 

E 

< 

? 

cc 

a> 

UT 

W 

a> 

> 

a> 

■o 

3 

*5. 

E 

< 

a> 

*o 

3 

*5. 

E 

< 

CO  £ 

TO 

w 

cn 

c 

o 

E 

TO 

TO 

L_ 

L» 

_ 

CSI  ^ 

co 

— 1 

L. 

ir\  i— 

TO 

nO  —i 

TO 

O-  — 1 

3 


E 

3 


15 

CONFIDENTIAL 


CONFIDENTIAL 


Assembly 


r/o* 


PfVSSLAF  PCLEF  POETS- 


CONFIDENTIAL 


Fig.  2  The  50-Inch  M«ch  8  Tunnel  (B) 
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Tunnel  Test  Section 
Fig.  3  The  50- Inch  Mach  10  Tunnel  (C) 
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Assembly 


Nozzle  and  Test  Section 
Fig.  4  The  12- Inch  Supersonic  Tunnel  (D) 
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Fig.  5  The  12-Inch  Hypersonic  Tunnel  (E) 
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Eig .  11  Bearing  Pivot,  Free  Oscillation,  Transverse  Rod  Supported 
Balance  System,  Mounted  in  the  50-Inch  Mach  8  Tunnel  (B) 


Eccentric 
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Fig.  12  Angular  Transducer  Bearing  Assembly 
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Mach  Number 

Fig.  14  Variation  of  the  Ratio  of  Model  Still-Air  Damping 
to  Model  Aerodynamic  Damping  with  Mach  Number 


Ref.  7 

12-Inch  Hypersonic 
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Fig.  15  Variation  of  the  Damping-in-Pitch  Derivatives  with  Mach  Number 
for  the  Pershing  Re-entry  Vehicle 
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A  SIMPLIFIED  INERTIA -COMPENSATED  BALANCE 
TECHNIQUE  FOR  WIND  TUNNEL  MEASUREMENT 
OF  LAUNCH  VEHICLE  RANDOM  BUFFET  EXCITATION 


by  C.  V.  Stahle,  C.  G.  Stouffer  and  W.  Silver 

Space  Systems  Division 
Martin  Company 
Baltimore  3,  Maryland 

Abstract 


A  simplified  technique  for  measurement  of  launch  vehicle  random  buffet 
excitation  is  presented  which  uses  a  rigid  model  of  the  forward  portion  of 
the  vehicle.  The  technique  combines  the  signals  from  a  standard  two-point 
moment  balance  with  the  signals  from  two  strategically  located  accelerometers 
to  provide  an  inertia  compensated  measurement  of  the  aerodynamic  moment 
at  two  points  on  the  model.  By  proper  design  of  the  model  and  balance,  the 
forward  accelerometer  will  completely  compensate  for  the  inertia  effects  at 
the  forward  balance  location,  and  the  aft  accelerometer  will  compensate 
for  the  inertia  effects  of  the  first  three  model-sting  modes  at  the  aft  balance 
location.  For  the  frequency  range  usually  investigated,  this  provides 
adequate  compensation  of  the  model  inertia  effects  permitting  the  moment 
of  oscillatory  aerodynamic  forces  to  be  measured  by  the  Inertia  Com¬ 
pensated  Balance. 

It  is  shown  that,  with  a  linearized  partial  mode  representation  of  the 
vehicle  bending  modes,  the  modal  excitation  can  be  approximated  by  the 
integrated  moment  of  oscillatory  aerodynamic  forces  over  the  forward 
portion  of  the  vehicle  about  the  nodal  points.  These  moments  are  determined 
from  the  balance  measurements.  Both  the  Power  Spectral  Density  of  the 
modal  excitation  and  the  Cross  Spectral  Density  of  the  excitation  between 
modes  are  obtained.  From  these  spectra,  the  mean  square  loads  can  be 
evaluated  using  an  extension  of  the  Thompson-Barton  approximation  which 
includes  the  effects  of  the  correlation  in  the  vehicle  response  resulting  from 
the  simultaneous  excitation  of  the  modes. 

The  adequacy  of  the  technique  is  verified  experimentally  by  comparing 
the  random  moment  applied  to  the  model  by  an  electromagnetic  shaker  with 
that  measured  by  the  balance.  The  close  comparison  of  NASA  Model  8  wind 
tunnel  buffet  measurements  with  previous  pressure  integration  measure¬ 
ments  verifies  the  adequacy  of  the  pressure  approach  over  the  frequency 
range  investigated.  The  Inertia  Compensated  Balance  technique  can  be  ex¬ 
tended  to  measure  simultaneously  the  excitation  in  two  planes  and  is  applic¬ 
able  to  other  similar  investigations. 
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Nomenclature 

Cross  spectral  density  of  the  aerodynamic  moment  about  the 
node  of  mode  n  and  mode  m 

Cross  spectral  density  of  the  generalized  modal  excitation  of 
mode  n  and  mode  m 

Structural  damping  coefficient  in  mode  n 

Constant  relating  the  inertia  moment  at  the  bridge  location  to 
the  acceleration  at  the  accelerometer  location 

Distance  from  the  nose  of  the  vehicle  to  the  node  of  mode  n 

Jm(x)  h^  (x)  dx,  the  equivalent  mass  for  the  nth  mode 

Moment  at  station  X  in  the  nth  mode 

Aerodynamic  moment  at  the  balance  location 

Inertia  moment  at  the  balance  location  due  to  model  vibration 

Total  moment  sensed  by  the  balance  bridge 

Fourier  transform  of  the  effective  pressure  force  acting  on  an 
axial  element  at  X 

Complex  conjugate  of  the  Fourier  transform  of  the  effective 
pressure  force  on  an  axial  element  at  X' 

Time,  used  to  denote  limiting  process 

Time,  used  to  denote  time  variations 

Power  spectral  density  of  the  bending  moment 

Power  spectral  density  of  the  generalized  modal  excitation 

Power  spectral  density  of  aerodynamic  moment  about  the  nodal 
point  of  mode  n 

Acceleration  of  the  model  at  the  accelerometer  location 
Distance  aft  of  the  nose  of  the  prototype  vehicle  or  the  model 
Impedance  of  the  nth  mode  of  the  prototype  vehicle 
Deflection  shape  of  the  nth  mode 

Effective  modal  slope  over  forward  portion  of  the  vehicle 
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,  2 

Tan  g  / 1  -  -tj—  the  phase  angle  associated  with  the  impedance 

u 

n 

of  the  nth  mode 

Resonant  frequency  of  the  nth  mode 
Modal  coordinate  indicating  time  variation 
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Introduction 


This  paper  presents  a  simplified  experimental  technique  for  the  measure¬ 
ment  of  the  random  buffet  excitation  of  launch  vehicle  transverse  bending 
modes  using  a  rigid  model.  As  a  launch  vehicle  passes  through  the  transonic 
region,  it  is  subjected  to  large  random  pressure  fluctuations  due  to  boundary 
layer  separation  and  shock  wave  interaction.  These  pressure  fluctuations 
occur  primarily  over  the  forward  portion  of  the  vehicle  and  particularly  at 
locations  of  abrupt  change  in  shape.  Because  the  effects  of  these  fluctuating 
pressures  cannot  be  predicted  to  the  required  degree  of  accuracy  by  ana¬ 
lytical  methods,  they  are  investigated  experimentally  by  wind  tunnel  tests 
using  scaled  models.  The  model  similitude  requirements  and  the  scaling  of 
the  test  results  have  been  the  subject  of  other  papers  presented  at  this  sym¬ 
posium  and  will  not  be  discussed  in  this  paper. 

There  are  two  fundamentally  different  approaches  to  the  wind  tunnel  in¬ 
vestigation  of  the  transonic  buffet  excitation  of  launch  vehicle  bending  modes. 
One  approach  is  to  use  dynamically  scaled  models;  the  other  is  to  measure 
separately  the  modal  excitation  due  to  fluctuating  pressures  and  the  motion 
dependent  unsteady  aerodynamic  effects.  This  paper  is  concerned  only  with 
the  first  part  of  the  latter  approach,  tne  measurement  of  the  modal  excita¬ 
tion  due  to  random  pressure  fluctuations. 

The  excitation  of  the  vehicle  bending  modes  is  defined  by  the  integral  of 
the  instantaneous  pressure  fluctuations  over  the  surface  of  the  vehicle, 
weighing  each  area  element  by  its  associated  modal  deflection.  The  Power 
Spectral  Density  and  Cross  Spectral  Density  of  these  integrated  pressure 
fluctuations  provide  the  desired  measure  of  the  modal  excitation.  The  ap¬ 
proach  used  in  previous  wind  tunnel  tests  has  been  to  determine  this  modal 
excitation  from  measurements  of  the  fluctuating  pressures  acting  on  the 
vehicle.  Although  pressure  fluctuation  measurements  are  necessary  for  the 
determination  of  local  structural  excitation,  these  measurements  are  not 
amenable  to  analysis  of  the  overall  vehicle  response.  It  becomes  an  im¬ 
possible  task  to  determine  the  modal  excitation  if  only  uncorrelated  pressure 
measurements  are  made. 

To  provide  more  meaningful  data  than  offered  by  the  measurement  of 
pressures  alone.  Cole  and  Coe  developed  a  pressure  integration  technique 
for  measurement  of  the  integrated  modal  excitation.  Ref.  1.  Their  tech¬ 
nique  utilized  the  signals  from  a  large  number  of  pressure  transducers  lo¬ 
cated  over  the  forward  portion  of  the  model.  By  assigning  effective  areas 
to  the  transducers  and  electronically  summing  these  signals,  the  integrated 
modal  excitation  was  obtained.  The  difficulties  with  this  technique  are  due 
to  the  approximations  made  in  assigning  effective  areas  to  the  transducers 
and  to  the  complexity  of  the  instrumentation  and  data  reduction  systems. 

A  large  number  of  pressure  transducers  and  associated  electronic  equip¬ 
ment  must  be  used  even  for  measurements  on  simple  bodies  of  revolution 
and  cannot  be  confidently  applied  to  more  complicated  payload  shapes  such 
as  the  Dyna-Soar.  Because  the  excitation  of  only  one  mode  at  a  time  can 
be  determined  with  this  system  without  doubling  the  already  complex  data 
reduction  system,  the  cross  spectra  of  the  excitation  between  the  various 
vehicle  modes  cannot  be  obtained. 
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Because  of  these  difficulties  and  the  necessary  approximations  made,  a 
simplified  technique  for  more  accurately  measuring  the  integrated  excitation 
of  the  vehicle  modes  and  the  cross  spectra  of  the  excitation  between  modes 
was  desired.  Extending  the  work  of  the  Cornell  Aeronautical  Laboratory 
(which  had  developed  an  Inertia-Compensated  Balance  for  measurement  of 
static  moments  in  a  hypersonic  shock  tunnel)  and  in  conjunction  with  the  Aero¬ 
space  Corporation,  an  Inertia -Compensated  Balance  (ICB)  technique  has  been 
developed  by  Martin  which  fulfills  these  requirements. 

The  present  Iner.tia-Compensated  Balance  technique  uses  the  combined 
signals  from  two  strategically  located  accelerometers  and  a  two -point  moment 
balance  to  measure  the  modal  excitation  of  the  vehicle  bending  modes.  The 
arrangement  of  the  model  and  the  instrumentation  is  shown  in  Fig.  1.  The 
accelerometers  provide  inertia  compensation.  Each  accelerometer  is  lo¬ 
cated  and  its  signal  sensitivity  is  adjusted  such  that  its  electrical  signal  is 
equivalent  to  the  inertia  bending  moment  at  its  associated  balance  bridge 
location.  Hence,  the  difference  between  the  signals  from  one  balance  bridge 
and  one  accelerometer  provides  a  measure  of  the  integrated  moment  of  the 
aerodynamic  pressures  acting  over  the  forward  portion  of  the  model,  com¬ 
pensated  for  the  model  inertia  effects.  From  the  simultaneous  measurement 
of  the  integrated  moment  of  the  fluctuating  pressures  about  two  points  on  the 
model,  the  Power  Spectral  Density  and  the  Cross  Spectral  Density  of  the 
modal  excitations  can  be  determined  using  a  linearized  partial  mode  repre¬ 
sentation  of  the  prototype  bending  modes  similar  to  that  of  Ref.  1. 

In  the  first  section  of  this  paper,  the  determination  of  the  modal  excita¬ 
tion  from  moment  measurements  will  be  discussed.  Then,  the  principle  of 
operation  of  the  Inertia-Compensated  Balance  and  the  final  adjustment  and 
calibration  of  the  ICB  system  at  the  wind  tunnel  will  be  presented.  Random 
vibration  test  results  verifying  the  adequacy  of  the  system  will  be  presented. 
Measurements  of  buffet  excitation  obtained  from  wind  tunnel  tests  using  the 
ICB  technique  will  be  compared  with  data  obtained  previously  by  NASA  Ames 
for  a  similar  model  using  the  pressure  integration  technique.  The  use  of  the 
ICB  measurements  for  the  determination  of  mean  square  loads  will  be  dis¬ 
cussed  and,  finally,  the  use  of  this  technique  for  similar  investigations  and 
possible  extensions  of  it  will  be  mentioned. 


BALANCE 

■PARTIAL 
MODE  MODEL 

FIC.  1.  TYPICAL  MODEL  ARRANGEMENT  FOR  INERTIA  COMPENSATED  BALANCE 
MEASUREMENTS 
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Determination  of  the  Modal  Excitations  from  Moment  Measurements 


The  vehicle  load  due  to  random  buffet  excitation  in  which  the  designer  is 
generally  interested  is  the  magnitude  of  the  mean  square  bending  moment. 

To  determine  this  load,  the  power  spectral  density  of  the  bending  moment  due 
to  the  response  of  the  vehicle  in  its  bending  modes  must  be  known.  This  bend¬ 
ing  moment  power  spectral  density,  W„w  (u,  x),  can  be  defined  by  the  follow¬ 
ing  equation  from  Ref.  2: 


W 


BM 


(u.  x)  =  y  y 


M'  (x)  M'  (x)  e 
n  m 


-i  (6 


0  ) 
m 


n  m 


^r 


f  vu* 


x)  c  *(u,  x1)  <j>  (x)  <)>  (x')  dx  dx1 

p  n  m 


(1) 


The  first  portion  of  the  right-hand  side  of  this  equation  represents  the  "bending 
moment  admittance"  of  the  vehicle's  structure  in  its  various  modes,  while  the 
second  portion  defines  the  generalized  modal  excitation.  The  modal  excitation 
due  to  random  forces  can  be  further  defined  as: 

WQ  l.u)  n  =  m 

C  C  if  F  (u,  x)  F  *(u,  x')  6  (x)  *  (x')  dx  dx'  =  n  (2) 

J  J  T  p  P  n  m  C  Q  (U)  n  4  m 

^n  m 

Hence,  the  Power  Spectral  Density  of  the  bending  moment  is  equal  to  the  sum 
over  all  modes  of  the  "benching  moment  admittance"  times  the  Power  Spectral 
Density  or  Cross  Spectral  Density  of  the  modal  excitation.  The  purpose  of  the 
ICB  technique  is  to  experimentally  evaluate  these  Spectral  Densities  of  the 
modal  excitation. 


The  integral  of  Eq  (2)  can  be  approximated  by  integrating  only  over  the  for¬ 
ward  portion  of  the  vehicle  because  the  pressure  fluctuations  and  the  modal 
deflections  are  generally  large  only  over  this  portion  of  the  vehicle.  In  addi¬ 
tion,  the  vehicle  mode  shape  can  be  represented  by  a  straight  line  for  the  for¬ 
ward  portion  of  the  vehicle  with  reasonable  accuracy  as  shown  in  Fig.  2.  With 
this  linearized  partial  mode  representation,  the  integral  of  Eq  (2)  can  be  written 
as; 


Fp(u'  x)  Fp*(u>  x>)  ^n 


4>' 

Tm 


X'  I  ,  ,  . 

j—  dx  dx'  = 

m  _ 


n 

d> ’  *'  C,,  (u) 

n  m 


n  =  m 
n  /  m 


(3) 
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Hence,  with  a  linearized  partial  mode  representation  of  the  full-scale  vehicle 
modes,  the  modal  excitation  is  proportional  to  the  Power  Spectral  Density  of 
the  moment  of  fluctuating  pressure  forces  about  the  nodal  points  of  the  modes 
and  the  Cross  Spectral  Density  of  the  moments  about  the  nodal  points  of  two 
different  modes. 


- ACTUAL  MODAL  D ERECTION 

- PARTIAL  MODE  REPRESENTATION 


£ 


Having  established  the  properties  of  the  buffet  pressure  fluctuations  that 
are  required  to  define  the  modal  excitation  using  a  linearized  partial  mode 
representation  of  the  prototype  bending  modes,  now  consider  the  capability 
of  the  ICB  system  to  measure  these  properties.  In  Appendix  A  of  this  paper, 
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it  is  shown  that  the  Power  Spectral  Density  and  Cross  Spectral  Density  of  the 
moments  at  the  two  ICB  balance  locations  can  be  used  to  determine  the  Power 
Spectral  Density  of  the  moment  about  any  other  point  as  well  as  the  Cross 
Spectral  Density  between  the  moments  at  any  two  points.  Because  the  modal 
excitations  have  been  defined  in  terms  of  the  Power  Spectral  Density  and  the 
Cross  Spectral  Density  of  the  moments  about  the  nodal  points  of  the  modes, 
these  modal  excitations  can  be  completely  determined  from  the  ICB  measure¬ 
ments  at  the  balance  locations. 

A  method  for  determining  the  mean  square  vehicle  loads  from  the  Power 
Spectral  Density  and  Cross  Spectral  Density  of  the  modal  excitations  is  pre¬ 
sented  in  Appendix  B  which  extends  the  Thompson -Barton  approximation  to 
include  the  load  contributions  due  to  the  simultaneous  excitation  of  the  vehicle 
modes. 


Principle  of  Operation  of  the  ICB  System 


If  a  standard  balance  without  inertia  compensation  were  used  to  measure 
the  moment  at  a  point  on  the  model,  the  balance  would  measure  both  the  in¬ 
tegrated  moment  of  pressure  fluctuations  and  the  moment  due  to  vibration  of 
the  model.  The  dynamic  sensitivity  of  the  balance  would  reflect  the  impedance 
of  the  model -sting  dynamic  system  having  peaks  at  model  resonances  and  low 
points  between  resonances  as  shown  in  Fig.  3.  If  the  model  inertia  effe  rts 
were  compensated,  the  calibration  curve  would  be  approximately  a  horizontal 
line  with  the  input  and  output  moments  of  the  balance  system  being  nearly  equal. 
This  type  of  compensation  is  accomplished  by  electronically  combining  an  ac¬ 
celerometer  signal  with  a  balance  signal  to  eliminate  the  inertia  effects  of 
model  vibration  on  the  measured  moment. 


FIG.  3.  DYNAMIC  SENSITIVITY  OF  UNCOMPENSATED  BALANCE 


The  principle  of  operation  of  the  ICB  system  can  be  explained  by  consider¬ 
ing  one  accelerometer -bridge  combination  of  the  typical  model  arrangement 
shown  in  Fig.  1.  The  moment  measured  by  the  strain  gage  bridge  of  the  bal- 
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ance  will  be  the  sum  of  the  aerodynamic  moment  and  the  moment  due  to  model 
inertia  effects.  This  can  be  written  as: 

MT(Xb,  t)  =  MA(Xb,  t)  +  Mj(Xb.  t)  (4) 


where  the  effects  of  damping  are  neglected  because  the  damping  is  generally 
small  and  would  affect  the  bridge  moment  only  in  narrow  frequency  bands 
near  the  model  resonances.  Considering  a  flexible  model,  the  inertia  moment 
at  the  bridge  location,  Xb,  can  be  written  in  terms  of  the  model  modes  and  the 


acceleration  at  a  point  on  the  model,  X  ,  as: 


W- 


w.(X  , 
1  a 


t) 


(5) 


where  w.(X.  ,  t)  is  the  modal  acceleration  at  X  . 

1  cl  cl 

If  by  proper  selection  of  an  accelerometer  location,  the  term  in  parentheses 
can  be  made  constant  for  all  modes,  Eq  (5)  can  be  written  as: 


Mj(Xb,  t)  =  K  )  w.(Xa,  t)  =  K  wT(Xa,  t)  (6) 


where 

MJ(Xb)  MJ(Xb*> 

K - 5 -  =  —7 -  =  ~2 -  =  .  .  .  (6a) 

"lW  “2W  "3*3<Xa> 

That  is,  the  inertia  moment  at  the  balance  bridge  will  be  proportional  to  the 
total  acceleration  at  a  point  on  the  model.  Then,  from  Eq  (4),  the  aerodynamic 
moment  can  be  obtained  by  subtracting  the  inertia  moment  which  is  propor¬ 
tional  to  the  model  acceleration  at  X  : 

a 

MA(Xb’  =  '  K  ^T(Xa'  t}  (7) 


Hence,  by  subtracting  the  acceleration  measured  by  a  properly  located  accel¬ 
erometer  and  multiplied  by  an  appropriate  "gain  constant,  "  K,  from  the  mo¬ 
ment  measured  by  the  balance,  the  moment  of  aerodynamic  pressure  fluctua¬ 
tions  is  obtained. 

The  instrumentation  system  used  to  accomplish  this  is  relatively  simple 
and  is  presented  in  Appendix  C.  Of  prime  importance  are  the  amplitude  and 
phase  characteristics  of  the  transducers  which  must  be  carefully  matched  to 
assure  adequate  compensation. 

If  only  one  mode  of  the  system  is  to  be  compensated,  the  requirement  of 
Eq  (6a)  can  be  relaxed  and  Eq  (6)  becomes: 
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MjIXj,  t) 


-  M;Kb)  - 

7i  W. 


Wj(Xa.  t)  =  Kj.  w1(Xa,  t) 


(6b) 


This  means  that  for  any  point  of  acceleration  measurement,  X  ,  the  inertia 

effects  due  to  model  motion  can  be  compensated  by  proper  determination  of  a 
"gain  constant,  "  K^,  relating  the  moment  at  the  balance  bridge  location,  X^, 

to  the  point  at  which  the  acceleration  is  measured.  This,  in  effect,  is  what 
was  accomplished  previously  by  the  Cornell  Aeronautical  Laboratory. 

In  order  to  provide  compensation  of  two  modes,  the  requirement  imposed 
is  that: 


M’i(xb>  .  “W 
"1  “2  *2<xa» 


(6c) 


For  a  given  model  and  balance  configuration,  this  relation  determines  the  lo¬ 
cation  of  the  accelerometer.  The  accelerometer  must  be  positioned  at  a  point 
such  that  the  ratio  of  modal  moment  at  the  balance  location  to  the  modal  ac¬ 
celeration  at  the  accelerometer  location  is  the  same  for  both  modes.  If  the 
modal  bending  moment  at  the  balance  bridge  location  is  considered  to  be 
fixed  (i.  e. ,  the  model  dynamic  characteristics  and  the  balance  geometry  are 
not  variable),  a  unique  accelerometer  location  is  determined. 


Extendin'*  the  compensation  to  include  three  modes,  the  requirement  in¬ 
dicated  by  L  (6a)  becomes 


M-(Xb)  M^(Xb)  M^(Xb) 

4  *l<Xa>  "  4  W  "  “3  W 


(6d) 


Now,  to  provide  compensation  of  three  modes,  the  modal  moment  at  the  bal¬ 
ance  location  must  be  a  variable.  During  the  design  of  the  model,  the  solu¬ 
tion  to  Eq  (6d)  can  be  determined  using  analytical  mode  shapes  to  evaluate  the 

modal  moment  at  the  bridge  location,  Mj(Xb).  The  accelerometer  and  bridge 

locations  at  which  the  ratio  of  modal  moment  to  modal  acceleration  are  equal 
for  the  first  and  second  modes  and  for  the  first  and  third  modes  are  deter¬ 
mined  for  the  model-balance  configuration  being  considered  using  relations 
similar  to  that  of  Eq  (6c).  The  results  are  plotted  as  shown  in  Fig.  4  and  the 
design  is  varied  until  the  first  three  modes  are  compensated  at  the  balance 
locations  as  indicated  in  the  figure. 

If  the  portion  of  the  model  forward  of  the  balance  bridge  location  is  effec¬ 
tively  rigid,  the  accelerometer  will  be  located  at  the  center -of -percussion  of 
the  model  about  the  bridge  location  and  the  inertia  effects  of  all  modes  will  be 
compensated.  By  proper  model  design,  this  can  be  accomplished  for  the  for¬ 
ward  balance  bridge  location  for  the  frequency  range  of  interest.  However, 
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the  flexibility  introduced  by  the  forward  balance  flexure  is  sufficient  to  in¬ 
validate  this  center-of -percussion  relation  for  the  aft  bridge.  Since  the  low 
frequency  range  is  of  interest,  compensation  of  only  the  first  three  model 
modes  will  generally  be  adequate  for  the  aft  balance  bridge.  This  is  ac¬ 
complished  by  varying  the  model  and  balance  configuration  as  discussed 
previously.  Although  guided  by  the  design  analysis  for  compensation  of  the 
first  three  model  modes,  final  adjustments  of  the  accelerometer  locations 
and  the  mass  of  the  model  are  made  at  the  wind  tunnel  to  achieve  optimum 
inertia  compensation. 


Final  Adjustment  and  Calibration 

After  assembling  the  model  and  balance  on  the  sting  in  the  wind  Jtunnel, 
final  adjustments  and  calibrations  of  the  Inertia  Compensated  Balance  system 
are  accomplished.  First,  the  accelerometers  are  located  for  a  nominal  mass 
of  the  model  and  the  adequacy  of  the  compensation  determined.  If  the  compen¬ 
sation  is  not  within  acceptable  limits,  the  model  mass  is  altered  by  the  addition 
or  removal  of  ballast  and  the  first  step  repeated.  Finally,  the  dynamic  sen¬ 
sitivity  of  the  balance  is  determined  o^er  the  frequency  range  of  interest.  The 
procedure  and  results  obtained  are  discussed  below. 

To  determine  the  optimum  accelerometer  location,  the  model-sting  system 
is  sinusoidally  excited  at  its  first  three  resonances  with  an  electromagnetic 
shaker  attached  to  the  sting  aft  of  the  balance.  With  the  model  excited  at  the 
sting  as  shown  in  Fig.  5,  the  moments  measured  by  the  balance  bridges  are 
due  entirely  to  model  inertia  forces.  With  the  system  excited  in  its  funda¬ 
mental  mode,  the  amplitudes  of  the  accelerometer  signals  are  adjusted  so  that 
the  difference  between  each  accelerometer  signal  and  its  corresponding  strain 
gage  signal  is  zero,  i.  e.  ,  the  amplitudes  of  the  accelerometer  signal  and  the 
bridge  signal  are  made  equal.  The  second  and  third  resonances  are  then  ex¬ 
cited  in  the  same  manner  and  the  signals  recorded.  The  accelerometers  are 
then  relocated  and  the  procedure  repeated.  The  optimum  locations  of  the  ac¬ 
celerometers  are  then  determined  by  plotting  the  ratio  of  acceleration  to  strain 
as  shown  in  Fig.  6.  As  mentioned  previously,  an  accelerometer  location  at 
which  this  ratio  is  equal  for  all  three  modes  is  desired.  Although  there  will 
be  some  response  due  to  nonresonant  modes  (i.  e. ,  first  mode  response  when 
the  system  is  excited  at  the  third  mode  resonance),  this  should  be  unimportant 
in  that  the  ratio  is  finally  made  equal  for  the  first  three  modes,  the  modes 
which  should  contribute  the  majority  of  nonresonant  response. 

If  the  compensation  is  within  acceptable  limits,  the  gains  and  locations  of 
the  accelerometers  are  fixed  and  remain  so  for  the  duration  of  the  test. 
Otherwise,  the  model  ballast  is  varied  and  the  preceding  steps  repeated  until 
adequate  compensation  is  obtained. 

With  the  proper  position  of  the  accelerometers  established  and  the  gains 
adjusted  and  locked  into  position,  the  dynamic  sensitivity  of  the  Inertia  Com¬ 
pensated  Balance  can  now  be  determined.  The  shaker  is  attached  to  the  model, 
forward  of  the  balance,  through  a  calibrated  force  transducer.  The  output  of 
the  ICB  system  (the  difference  between  the  balance  bridge  and  accelerometer 
signals)  should  measure  only  the  oscillatory  moment  applied  to  the  model, 
since  inertia  effects  should  be  cancelled.  Hence,  by  sinusoidally  forcing  the 
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model  through  the  frequency  range  of  interest,  the  dynamic  sensitivity  of  the 
system  can  be  obtained. 

Typical  results  are  shown  in  Fig.  7,  which  indicates  that  the  sensitivity  of 
the  system  is  constant  well  above  the  first  three  resonant  modes  of  the  model 
and  sting  support  system.  A  maximum  deviation  of  10%  is  indicated  through¬ 
out  the  frequency  range  of  10  to  140  cps,  which  was  the  range  of  interest  for 
this  particular  model.  Similar  results  were  obtained  for  both  sums  for  two 
different  shaker  locations.  In  contrast  to  the  ICB  system,  the  inability  of  an 
uncompensated  balance  to  measure  the  sinusoidal  input  moment  was  shown  in 
Fig.  3. 


Verification  of  the  ICB  System 

A  final  verification  of  the  adequacy  of  the  ICB  system  to  compensate  for 
inertia  effects  was  accomplished  by  applying  a  random  excitation  force  to  the 
model  forward  of  the  balance  through  the  calibrated  force  transducer.  A  typ¬ 
ical  comparison  of  the  Power  Spectral  Density  of  the  random  input  moment 
measured  with  the  force  transducer  and  the  output  moment  measured  with  the 
Inertia  Compensated  Balance  is  shown  in  Fig.  8.  Similar  results  were  ob¬ 
tained  for  both  balance  locations  with  excitation  at  two  different  model  loca¬ 
tions.  The  results  indicated  the  capability  of  the  ICB  system  to  measure  the 
Power  Spectral  Density  of  the  applied  random  moment  within  approximately 
±  1  db.  The  adequacy  of  the  system  calibration  is  also  verified  by  this  com¬ 
parison. 


Comparison  with  Pressure  Integration  Measurements 

Having  proven  the  capability  of  the  ICB  system  to  compensate  for  model 
inertia  effects,  wind  tunnel  tests  were  conducted  to  obtain  buffet  excitation 
measurements  comparable  to  those  obtained  previously  by  NASA  Ames  using 
the  pressure  integration  technique.  For  this  comparison,  the  NASA  Model  8 
hammerhead  configuration  was  selected  because  of  the  large  buffet  excitation 
measured  by  Ames.  The  model  was  of  the  same  scale  and  was  tested  in  the 
same  wind  tunnel  facility,  the  Ames  14 -foot  transonic  tunnel,  as  was  the  modei 
used  for  the  pressure  integrc'-ion  measurements. 

The  model  configurations  are  shown  in  Fig.  9,  which  compares  the  two 
measurement  systems.  The  pressure  integration  measurements  utilized  the 
signals  from  76  pressure  transducers  located  at  17  stations  over  the  forward 
portion  of  the  model.  Four  or  eight  pressure  transducers  were  used  at  each 
station  and  combined  in  a  bridge  to  give  the  moment  contribution  over  a  small 
span.  The  transducer  locations  were  selected  to  provide  approximately  equal 
area  moment  contributions  about  station  64.  85,  one  point  about  which  the  mo¬ 
ment  was  measured.  The  signals  from  each  section  were  then  combined 
through  a  series  of  operational  amplifiers  to  give  the  total  integrated  time  cor¬ 
related  moment  about  either  of  the  two  reference  points  indicated  in  Fig.  9. 

The  frequency  range  of  interest  was  from  approximately  20  to  140  cps,  which 
would  cover  the  simulated  frequency  range  of  structural  bending  modes  for 
typical  flight  vehicles. 
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FIG.  6.  SURVEY  OF  ACCELEROMETER  LOCATIONS  FOR  BOTH  MODEL  STRAINGAGES 


v>l 


na<* 

o  />  O 

0 

LU 

a 

1 

FORWARD  BRIDGE 

S  i 

a 

Q£ 

rr, 

1  1. 

I  , 

, 

, 

0  20 

40 

60  80  100 

120 

140 

t—  | 

3l 


1.2 
1.0 
0.8 
16 
0. 4i 
12 


Uh°^ 


g  §  § 

1-Sis 
J _ L _ L 


AFT  BRIDGE 


20 


40  60  80  100 

FREQUENCY  (CPS) 


120  140 


FIG.  7.  DYNAMIC  SENSITIVITY  OF  INERTIA  COMPENSATED  BALANCE 
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FIG.  8.  COMPARISON  OF  RANDOM  MOMENT  PSD  MEASURED  WITH  THE 
INERTIA  COMPENSATED  BALANCE  AND  A  FORCE  TRANSDUCER 


The  relative  simplicity  of  the  Inertia  Compensated  Balance  is  evident  from 
Fig.  9.  Only  four  transducers  are  required  to  measure  simultaneously  the 
moment  at  two  model  locations  as  compared  to  a  single  location  with  the  pres¬ 
sure  integration  technique.  In  order  to  be  assured  that  the  buffet  excitation 
was  the  same,  pressure  measurements  were  made  at  the  three  model  locations 
indicated  in  Fig.  9  for  comparison  with  previous  measurements. 

The  Power  Spectral  Density  of  the  aerodynamic  moment  about  the  forward 
reference  point.  Station  43.  25,  is  shown  in  Fig.  10  for  Mach  0.  90.  Comparison 
of  the  ICB  data  and  the  pressure  integration  data  indicates  negligible  differences 
in  the  low  frequency  range  (10  to  80  cps)  and  a  maximum  difference  of  2  db  at 
110  cps.  This  is  considered  to  be  well  within  the  expected  accuracies  of  the 
data  considering  the  necessary  approximations  included  in  the  pressure  integra¬ 
tion  technique  and  could  possibly  be  due  to  differences  in  the  pressure  fluctua¬ 
tions  as  indicated  by  the  comparison  of  pressure  measurements  shown  in  Fig. 

11.  The  Power  Spectral  Density  of  the  pressure  fluctuations  at  this  measure¬ 
ment  location  also  differed  by  approximately  2  db  at  110  cps.  Good  repeatibil- 
ity  of  the  ICB  data  is  evident  in  Fig.  10  which  includes  data  points  obtained 
from  two  test  runs.  Similar  comparisons  were  also  made  at  Mach  0.  95  with 
comparable  results. 

A  comparison  of  the  Power  Spectral  Density  of  the  aerodynamic  moment 
at  the  aft  end  of  the  model,  Station  64.  85,  is  shown  in  Fig.  12  and  was  ob- 
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FIG.  9.  COMPARISON  OF  ICB  INSTRUMENTATION  AND  PRESSURE  INTEGRATION  INSTRUMENTATION 
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FIG.  11.  COMPARISON  OF  PRESSURE  POWER  SPECTRAL  DENSITY  AT 
STATION  32.60,  MACH  0.90 
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FIG.  12  COMPARISON  OF  AERODYNAMIC  MOMENT  PSD  AT  STATION 
64. 85,  MACH  0  90 
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tained  by  transferring  the  ICB  data  to  that  point.  Although  the  comparison 
with  the  pressure  integration  data  is  not  as  close  as  for  Station  43.  25  (having 
larger  differences  below  30  cps  and  above  120  cps),  it  is  considered  to  com¬ 
pare  favorably. 

On  the  basis  of  these  tests,  it  is  evident  that  the  Inertia  Compensated  Bal¬ 
ance  technique  provides  a  considerable  reduction  in  the  complexity  of  the  in¬ 
strumentation  even  for  a  simple  body  of  revolution.  The  comparison  of  the 
buffet  excitation  as  measured  by  the  two  techniques  provides  increased  con¬ 
fidence  in  the  ICB  system  and  verifies  the  capability  of  the  pressure  integra¬ 
tion  technique  to  measure  the  integrated  buffet  excitation  over  the  frequency 
range  investigated.  For  more  complicated  shapes  such  as  a  glider  configura¬ 
tion,  the  Inertia  Compensated  Balance  provides  the  only  practical  approach. 

In  addition,  more  meaningful  data  is  obtained  with  the  Inertia  Compensated 
Balance.since  the  simultaneous  measurement  of  the  moment  at  two  points  is 
obtained  enabling  the  modal  excitations  to  be  completely  determined  as  dis¬ 
cussed  in  the  first  section  of  this  paper. 


Conclusions 


On  the  basis  of  the  theoretical  and  experimental  results  of  this  investiga¬ 
tion,  the  following  conclusions  can  be  made  regarding  the  Inertia  Compensated 
Balance  technique: 

(1)  A  relatively  simple  experimental  technique  has  been  developed  for  the 
wind  tunnel  measurement  of  the  modal  excitations  due  to  random  buffet  pres  - 
sures  which  reduces  by  more  than  an  order  of  magnitude  the  number  of  trans¬ 
ducers  and  amplifiers  required. 

(2)  This  technique  enables  the  simultaneous  determination  of  the  excita¬ 
tion  of  several  modes  as  well  as  the  correlation  of  the  excitation  between 
modes. 

(3)  The  adequacy  of  this  technique  has  been  verified  by  several  vibration 
tests  of  a  model  indicating  an  accuracy  of  approximately  ±  1  db  for  the  Power 
Spectral  Density  of  the  measured  moment. 

(4)  Wind  tunnel  measurements  of  the  NASA  Model  8  buffet  excitation  agree 
almost  identically  with  previous  measurements  using  the  pressure  integration 
technique  verifying  the  adequacy  of  that  technique  over  the  frequency  range  in¬ 
vestigated. 

(5)  A  method  of  analysis  has  been  indicated  for  use  with  the  experimental 
data  obtained  and  includes  the  correlation  effects  between  modes  in  one  plane. 


Extensions  and  Applications 


This  paper  has  been  directed  at  the  specific  problem  of  determining  launch 
vehicle  transverse  bending  loads  due  to  transonic  buffeting  and  has  considered 
loads  in  only  one  plane.  However,  the  ICB  technique  can  be  applied  to  other 
investigations  which  impose  similar  measurement  requirements  and  can 
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readily  be  extended  to  include  two  planes  of  measurement.  Generally,  this 
technique  is  applicable  to  any  investigation  in  which  the  integrated  effects  of 
random  excitation  are  to  be  determined  for  modes  that  can  be  linearized. 
Wind  tunnel  investigations  of  oscillatory  ground  wind  loads  present  an  almost 
identical  measurement  problem  and  could  be  investigated  with  this  technique. 
By  simply  duplicating  the  balance  system  in  an  orthogonal  plane,  the  modal 
excitations  in  hoth  the  pitch  and  yaw  planes  could  be  investigated  simultane¬ 
ously.  This  would  enable  the  correlation  of  the  modal  excitations  between 
planes  to  be  determined,  which  would  be  of  prime  importance  for  vehicles 
having  axial  symmetry  whe  the  *  iodal  frequencies  in  two  directions  are 
equal.  Hence,  the  Inertia  Compensated  Balance  technique  can  be  used  for 
investigations  other  than  buffet  and  can  be  extended  to  measure  the  correla¬ 
tion  between  modal  excitations  in  two  directions. 
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Appendix  A 


Determination  of  the  Power  Spectral  Density  and  Cross  Spectral, Density 
of  the  Modal  Excitation  from  ICB  Measurements 

The  ICB  system  provides  the  moment  at  two  locations  as  indicated  in  Fig. 
A-l.  ; 


The  modal  excitation  has  been  defined  in  terms  of  the  Power  Spectral  Den¬ 
sity  and  Cross  Spectral  Density  of  the  moment  about  the  prototype  nodal 
points.  Because  the  points  at  which  the  moments  are  measured  with  the  ICB 
system  do  not  necessarily  correspond  to  the  prototype  nodal  points,  it  is  nec¬ 
essary  to  define  the  spectral  densities  at  the  nodal  points  in  terms  of  the 
spectral  densities  at  the  measurement  locations.  This  is  developed  in  this 
appendix  for  two  arbitrary  nodal  points,  X^  and  Xg. 

The  moment  at  any  point,  X,  can  be  defined  in  terms  of  the  measured  mo¬ 
ments,  Mp(t)  and  M^(t)  as: 

MX(t)  =  Gj+  X]  MR(t)  -  [  I  ]  MF(t)  <A'X) 

where  X  is  the  distance  aft  of  the  rear  measurement  point  as  indicated  in  Fig. 
A-l. 


The  Power  Spectral  Density  and  Cross  Spectral  Density  for  the  two  meas^ 


urement  points  are 

defined  as 

- 

Wp(f)  =  ^ 

FT*(f)  FT(f) 

(At2a) 

T 

WR(f)  =  t“„ 

RT*(f)  RT(f) 

(A -2b) 

T 

CFR<f>  =  TZ 

FT*(f)  RT(f) 

(A-2c) 

T 
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where 

T 

FT<f)  =  C 
-T 

T 

RT(f)  »  J 

and  the  asterisk  denotes  the  complex  conjugate. 

Then  the  Power  Spectral  Density  at  any  arbitrary  nodal  point,  X,  can  be 
written  as: 

WX(f>  ■  Tl"  J£  +  l]  RT*(t>  -  [J]  FT*«>(  j[‘  +  l]  RT(f> 

-  [?]  FT(t>|  <A'3i> 

'  [‘  +  l]  2  WR(f)  +  [|]  2  WF(fl  '  ['  +  a]  B]  CFR(f) 

-  [1  +  l]  [|]  CRF(t>  (A-3b» 

but, 

CFR(f)  ■  CRF*<» 

therefore, 

wx'»  ■  [l  +  a]  WR<»  +  [|]  wF<f)  - 2  [■  +  a]  [|]  cFR(f:  real> 

(A-3c) 

Consequently,  the  Power  Spectral  Density  at  any  arbitrary  nodal  point  can  be 
determined  from  the  Power  Spectral  Densities  of  the  measured  moments  and 
the  real  portion  of  the  Cross  Spectral  Density  (Co-Spectrum)  of  the  measured 
moments. 

Similarly,  the  Cross  Spectral  Density  between  any  two  nodal  points,  X^ 
and  Xg,  can  be  determined  in  terms  of  the  measured  spectra  by  defining 


At(I)  -  [l  +^l]  RT<f>  -  [^p]  Ft«> 

(A-4a) 

%•«>-  [>  *t]  RT<»  -  hrK<« 

(A- 4b) 
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Then,  the  Cross  Spectral  Density  between  the  moments  at  X^,  and  Xg  be¬ 
comes,  by  definition:  “  16 


CY  Y  (f)  = 
X1X2 


lim  AT*W»T«> 
T-m  T 


(A- 5a) 


lim 
T 


It  ♦*] 


RT(f) 


(A-5b) 


£,  T  |[1  +  t]rt*(i>  -  hr|v(tl 

-  m  rt<«i 

[4]  H 

-  [‘  +  ^r]  hr]  cRF<n  -  [l  +  t]  &]  cfr<«  <A-5c» 


WF<f> 


Separating  the  Cross  Spectral  Density  into  its  real  and  imaginary  parts  (the 
Co-Spectra  and  Quad-Spectra),  this  can  be  written  as: 


wF(f) 

real)  (A-5d) 

imaginary) 


From  the  relations  developed  in  Eqs  (A-3c)  and  (A-5d),  it  can  be  seen  that 
the  Power  Spectral  Density  at  any  nodal  point,  X,  and  the  Cross  Spectral  Den¬ 
sity  between  two  arbitrarily  selected  nodal  points,  Xj  and  Xg,  can  be  deter¬ 
mined  from  the  moment  Power  Spectral  Densities  at  the  balance  bridge  loca¬ 
tions  and  the  Cross  Spectral  Density  of  the  measured  balance  moments. 


Appendix  B 


Determination  of  the  Mean  Square  Loads 

Having  established  the  capability  of  the  ICB  system  to  measure  the  Power 
Spectral  Density  and  Cross  Spectral  Density  of  the  modal  excitation,  the  de¬ 
termination  of  the  mean  square  vehicle  loads  from  these  data  will  be  discussed. 
The  mean  square  load  at  each  vehicle  location  is  determined  from  the  integral, 
over  frequency,  of  the  Power  Spectral  Density  of  the  load  of  interest.  Because 
the  bending  moment  is  usually  the  desired  load,  the  following  discussion  will 
be  directed  toward  the  determination  of  the  mean  square  bending  moment  al¬ 
though  the  same  approach  can  be  used  for  other  loads.  By  combining  Eqs  (1) 
and  (2),  the  bending  moment  Power  Spectral  Density  can  be  written  as: 
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(B-l) 


In  this  equation,  the  Power  Spectral  Density  and  Cross  Spectral  Density  of 
the  excitation  are  determined  from  the  Inertia  Compensated  Balance  meas¬ 
urements  using  Eq  (3),  The  modal  moment  coefficients  are  determined  from 
the  vibration  analysis  of  the  vehicle,  and  the  modal  impedance  is  determined 
from  the  vibration  analysis  including  the  effects  of  motion -dependent  aero¬ 
dynamic  forces  (i,  e.,  aerodynamic  damping,  mass,  and  stiffness  effects). 

The  mean  square  bending  moment  is  determined  by  the  frequency  integral 
of  the  bending  moment  Power  Spectral  Density  as: 


M 


£  [m;<x|]  2  J 


w  u 

n 


Zn(u>) 


du 


II 

n  m 
n  /  m 


M»  (x)  M'  (x) 
n  m 
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-i  (8  -  8  ) 
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MTP 
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TST]  CQnQm(lj)  dw 


(B-2) 


For  small  damping  and  significantly  different  modal  frequencies,  this  inte¬ 
gral  can  be  approximated  as 

K(x)] 2  wq  s> 


M2(x) 


■51 


m"2  3  ‘ 
g  M  u 
6n  n  n 


»  V  V  Mn(x)  Mmlx) 

Z  M  M  u  <u2  -  u2) 

n  m  n  m  n  m  n' 


[cQ  Q  (un:  imaginary> 
L  v*ntm 


(J  I 

^  imaginary)l 


J 


(B-3) 


The  first  term  of  this  equation  is  the  Thompson-Barton  approximation.  Ref.  3, 
reflecting  the  resonant  response  in  the  various  vehicle  modes.  The  second 
term  reflects  the  correlation  between  the  response  in  two  modes  due  to  the  cor- 
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relation  of  the  modal  excitations  and  can  be  derived  in  a  similar  manner  as 
the  first  term. 


Appendix  C 


Inertia  Compensated  Balance  Instrumentation 


The  instrumentation  used  in  the  ICB  system  is  depicted  in  Fig.  C-l. 


TO 

TAPE 

RECORDER 


FIG.  C-l.  TYPICAL  ICB  CIRCUIT  SCHEMATIC 


The  symbols  used  indicate  the  following: 

Aj  Forward  accelerometer  signal 
Ar  Aft  accelerometer  signal 

Forward  strain  gage  bridge  signal 
Sr  Aft  strain  gage  bridge  signal 

Ej.  ICB  moment  measured  at  the  forward  bridge  location 
Er  ICB  moment  measured  at  the  aft  bridge  location 

Of  these.  A.,  S-,  A  ,  and  S  are  primary  measurements;  E,  and  E  are 
i  i  r  r  i  r 

produced  by  summations  of  their  respective  accelerometer  and  strain  gage 

bridge  signals. 

Sj  and  were  obtained  from  a  standard  two-point  balance  instrumented 
with  four  active  arm  bending  bridges.  The  strain  gages  used  were  Baldwin 
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Lima  Hamilton  Type  AB7.  These  bridges  were  excited  and  the  signals  were 
amplified  by  a  CEC,  Type  1-127,  20-kc  carrier  amplifier. 

Aj  and  Ar  were  obtained  using  Endevco  Type  2235  piezo-electric  accel¬ 
erometers,  Signals  from  these  transducers  were  amplified  and  matched  to 
the  remaining  system  by  Endevco  Type  2617  amplifiers. 

The  choice  of  a  piezo-electric  accelerometer  was  dictated  by  the  need  for 
a  transducer  with  a  high  resonance  frequency  and  low  damping.  This  is  re¬ 
quired  to  avoid  a  large  phase  shift  at  high  frequencies  between  the  actual  and 
indicated  acceleration.  At  low  frequencies,  a  phase  difference  existed  due  to 
the  roll-off  characteristic  of  the  accelerometer-amplifier  combination  al¬ 
though  no  amplitude  roll-off  existed  in  the  frequency  range  of  interest.  This 
phase  difference  was  removed  by  determining  the  roll-off  and  matching  it  in 
the  bridge  circuit  with  a  high  pass  filter.  With  this  device,  the  phase  differ¬ 
ences  in  the  strain  and  accelerometer  circuits  were  eliminated  without  causing 
an  amplitude  roll-off  by  adjusting  the  filter  corner  frequency  well  below  the 
frequency  range  being  investigated.  In  addition,  Af  and  accelerometer- 

amplifier  combinations  were  matched  in  phase  roll-off  characteristics  so  that 
instrumentation  phase  shifts  between  Ef  and  Ef  were  kept  negligibly  small. 

The  amplitude  effects  of  roll-off  were  outside  the  frequency  range  of  interest 
for  ail  components. 

Amplitude  adjustments  were  made  by  fixing  the  and  S  signals  at  con¬ 
venient  levels  and  adjusting  the  amplitudes  of  A  and  A  as  required  during 
the  calibration  of  the  system.  f  r 

Matched  low  pass  filters  were  used  in  the  Af,  Ar,  S^,  and  Sr  circuits. 

These  eliminated  undesired  high  frequency  signals  from  the  data.  These  sig¬ 
nals  were  primarily  20-kc  noise  from  the  carrier  amplifiers  and  acoustically 
induced  accelerations  during  the  wind  tunnel  tests.  Care  was  taken  to  assure 
that  these  filters  did  not  produce  undesirable  phase  and  amplitude  effects. 

Figure  C-l  shows  several  operational  ampli^rs  other  than  the  summing 
amplifiers.  These  were  used  primarily  as  impedance  matching  devices  for 
the  circuits  and  also  as  points  where  the  various  signals  could  be  monitored 
during  the  calibrations.  Both  the  operational  an.d  summing  amplifiers  were 
fabricated  by  NASA-Ames  personnel  from  commercial  components.  The  ac¬ 
celerometer  and  strain  gage  bridge  were  wired  so  that  the  strain  signals 
were  180°  out  of  phase  with  the  acceleration  signals.  These  signals  were 
then  summed. 

All  data  was  recorded  on  a  Precision  Instrument  Type  PI  214  14-track 
FM  Type  Recorder-Reproducer.  Calibrations  were  made  to  determine  the 
phase  shifts  inherent  between  critical  tracks  of  recorded  tape.  These  tended 
to  increase  with  frequency  and  were  held  to  less  than  2°  in  the  range  of  in¬ 
terest.  Two -minute  records  were  taken  at  each  data  point. 


Data  reduction  was  performed  by  the  Dynamics  -Laboratory  at  Martin- 
Baltimore.  The  data  were  reproduced  on  an  Ampex  Type  FL2Q0  Tape  Loop 
Reproducer.  Power  Spectral  Density  and  Cross  Spectral  Density  plots  were 
produced  using  a  Technical  Products  Company  Type  627  Analyzer  with  6  45 
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THE  EFFECT  OF  MODEL  SCALE  ON  RIGID -BOOT  UNSTEADY 
PRESSURES  ASSOCIATED  WITH  BUFFETING 
By  Charles  F.  Coe 

National  Aeronautics  and  Space  Administration 
Ames  Research  Center 
Moffett  Field,  Calif. 


ABSTRACT 


The  question  of  the  effect  of  model  scale  on  launch  vehicle  dynamic 
measurements  is  one  which  invariably  arises  in  connection  with  the  applica¬ 
tion  of  measured  unsteady  pressures  on  wind-tunnel  models.  The  scaling  of 
unsteady  pressure  measurements  is  discussed  in  this  paper  which  presents 
comparisons  of  results  of  pressure -fluctuation  measurements  on  both  wind- 
tunnel  models  and  on  full-scale  Ranger  5  and  Mercury  vehicles.  In  addition, 
results  of  tests  using  different  sized  models  are  shown  in  order  to  cover 
some  of  the  different  types  of  local  flow  associated  with  buffeting.  The 
effects  of  scale  on  the  root -mean -square  fluctuations  of  pressure,  the 
longitudinal  correlation  of  the  fluctuations,  and  on  the  power  spectral 
densities  are  shown  for  selected  transonic  Mach  numbers  where  the  fluctua¬ 
tions  are  most  severe. 
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THE  EFFECT  OF  MODEL  SCALE  ON  RIGID-BODY  UNSTEADY 
HUESUKES  ASSOCIATED  WITH  BUFFETING 

By  Charles  F.  Coe 

National  Aeronautics  and  Space  Administration 
Ames  Besearch  Center 
Moffett  Field,  Calif. 


INTRODUCTION 


The  unsteady  aerodynamic  loads  on  launch  vehicles  at  transonic  speeds 
have  recently  received  considerable  attention.  These  unsteady  loads  result 
from  pressure  fluctuations  that  occur  within  regions  of  shock  waves  or 
regions  of  separated  flow  or  both.  Escape  rocket  systems,  blunt  noses  to 
insure  abort  stability,  bulbous  payloads  which  are  larger  than  their  booster 
rockets,  and  numerous  protuberances  all  lead  to  this  troublesome  flow 
unsteadiness  which  can  cause  buffeting. 

Several  wind-tunnel  investigations  (refs.  1  to  6)  have  been  undertaken 
at  Ames  and  Langley  Research  Centers  and  at  Arnold  Engineering  Development 
Center  to  measure  pressure  fluctuations  cm  both  specific  corrf igurations  and 
also  on  a  variety  of  body  shapes  to  determine  effects  of  profile  variations. 

When  any  of  these  pres  sure -fluctuation  measurements  are  used  for  estimation 
of  vehicle  dynamic  response,  questions  invariably  arise  as  to  the  proper 
method  of  scaling  the  data  to  full  scale.  For  example,  are  the  concepts  of 
a  constant  Stroubal  number  and  scaling  by  application  of  the  commonly  used 
reduced  frequency  parameter,  wd/V,  appropriate  for  random  nonperiodic  buffet 
pressures.  For  lack  of  any  verifying  Information  the  above  concepts  have 
been  enqployed  in  references  6  and  7.  'Two  recent  investigations  (refs.  8 
and  9)  have  devoted  some  attention  to  the  problem  of  scaling  unsteady  pres¬ 
sures.  The  latter  (ref.  9),  which  contains  data  for  models  varying  in  size 
by  a  ratio  of  5  to  1,  tends  to  substantiate  the  validity  of  the  reduced  fre¬ 
quency  parameter  for  scaling.*  The  ultimate  test  of  scaling  buffet  pressures, 
however,  comes  with  the  direct  comparison  between  wind-tunnel  model  data  and 
full-scale  data  obtained  during  the  launching  of  a  vehicle. 

The  measurement  of  differential  pressure  fluctuations  at  two  Agena 
stations  during  the  Ranger  5  launch  and  of  pressure  fluctuations  on  the  Mercury- 
Atlas  adapter  during  launch  of  the  MA-4,  MA-5,  MA-7,  and  KA-S  provided  a  good 
opportunity  to  make  such  comparisons.  As  a  result,  tests  of  a  7 -percent  scale 
Mercury  model  and  a  10-percent  scale  Ranger  model  were  included  in  an  investi¬ 
gation  of  scale  effects  using  different  sized  models  which  was  in  progress  at 
Ames  Research  Center.  The  tests  were  conducted  at  transonic  Mach  numbers  with 


*Results  from  this  investigation  appear  in  the  paper  by  Hanson  and  Jones 
of  Langley  presented  at  an  earlier  session  of  this  Symposium. 
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pressure  transducers  located  at  the  same  stations  as  the  full-scale  flight 
instrumentation.  It  is  primarily  the  results  of  these  comparative  tests 
that  are  contained  herein.  Additional  data  from  the  general  research  pro¬ 
gram  are  included  to  illustrate  the  effects  of  scale  for  different  types 
of  unsteady  flow. 
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MODELS 


The  models  tested  to  Investigate  scaling  effects  (fig.  l)  vere  a 
7 -percent  scale  Mercury -Atlas,  a  10-percent  scale  Ranger,  and  three  sizes 
each  of  models  8  and  13a.*  For  the  Mercury -Atlas  and  Ranger  models  the 
Atlas  cylindrical  hodies  were  neither  scaled  in  length  nor  did  they  include 
any  protuberances.  Key  longitudinal  stations  are  indicated  in  the  figure  to 
aid  in  visualizing  the  placement  of  the  pressure  transducers  relative  to  the 
tody  profiles.  The  reference  diameters  used  for  scaling  the  results  are 
also  shown. 

All  the  models  were  sting  mounted  except  the  26-inch-diameter  model  13a 
which  was  tested  as  a  half  model  mounted  on  the  ll-foot  transonic  wind  tunnel 
wall  to  minimize  the  effects  of  blockage.  All  the  models  were  tested  in  the 
14-foot  wind  tunnel  except  the  Mercury -Atlas  which  was  tested  in  the  11 -foot 
transonic  wind  tunnel  and  the  6. 31 -inch-diameter  model  8  which  was  tested  In 
the  6-  by  6-foot  wind  tunnel. 

The  pressure  transducer  stations  that  match  the  locations  on  the  full- 
scale  vehicles  are  marked  with  an  x  in  figure  1.  The  differential  pressure 
fluctuations  were  measured  on  the  Ranger  by  taking  the  electrical  difference 
between  the  outputs  from  the  transducers  mounted  opposite  each  other.  Flush- 
mounted  l/4-inch-diameter  strain -gage -type  transducers  were  used  on  all  the 
models . 

A  photograph  of  the  Mercury  model  in  the  11-  by  11 -foot  wind  tunnel  is 
shown  in  figure  2.  An  arrow  indicates  the  relative  location  of  the  trans¬ 
ducer  with  respect  to  the  bulges  on  the  ring  clamp  between  the  spacecraft  and 
adapter  section.  The  transducer  station  behind  one  of  the  bulges  at  approxi¬ 
mately  mid-length  of  the  adapter  is  that  at  which  flight  data  were  obtained. 


EFFECTS  OF  SCALE 


When  dynamic  data  associated  with  motion  effects  are  measured  on  models, 
it  is  accepted  that  the  geometric  similarity  of  the  flow  cannot  be  maintained 
unless  the  reduced  frequencies  for  model  and  full  scale  are  the  same.  Although 
there  has  been  doubt  whether  the  reduced  frequency  parameter  is  applicable  to 
random  unsteady  aerodynamic  measurements,  simple  dimensional  analyses  by 
Liepmann  (ref.  10)  of  the  variable  parameters  involved  indicate  that  the  dimen¬ 
sionless  reduced  frequency  parameter  can  be  used  to  scale  pres sure -fluctuation 
measurements : 


md^  Cixdg 

vT =  v7 


*The  model  numbers  were  assigned  to  the  series  of  launch  vehicle  payload 
shapes  tested  at  Ames  Research  Center. 
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it  follows  that  in  order  to  compare  RMS  measurements  from,  model  to  full  scale 
or  from  one  model  size  to  another,  the  range  of  frequencies  included  in  the 
RMS  measurements  must  he  scaled  in  proportion  to  the  reference  dimensions. 


Comparisons  of  Mercury  Results 


A  comparison  of  power  spectra  of  the  pressure  fluctuations  measured  on 
the  Mercury -Atlas  adapter  is  shown  in  figure  3-  The  wind-tunnel  data  are  for 
a  fixed  Mach  number  of  1.0  while  the  flight  data  were  obtained  for  a  range  of 
Mach  numbers  near  1.*  The  comparison  is  based  on  reduced  frequency .  (The 

V/2itd _ in  the  ordinate  scale  converts  the  power  spectra  from  £Cp2  per  cps 

to  £Cp2  per  unit  of  reduced  frequency. )  It  can  be  seen  from  these  results 
that  the  7 -percent  model  data,  which  have  been  scaled  by  a  factor  of  about 
14. 3,  fit  reasonably  well  within  the  limits  of  the  spread  of  flight  measure¬ 
ments.  To  illustrate  the  extent  of  frequency  scaling  the  limit  of  the 
7 -percent  model  data  is  at  approximately  5,200  cps  while  the  full-scale  data 
for  the  same  reduced  frequency  is  approximately  365  cps. 

A  curve  is  also  shown  from  reference  6.  While  this  curve  appears  higher 
than  the  others,  the  difference  between  it  and  MA-4  levels  is  less  than  the 
spread  from  MA-4  to  MA-5-  This  spread  of  a  factor  of  about  4  on  a  mean-square 
scale  (factor  of  2  on  a  RMS  scale)  serves  to  indicate  the  limits  of  accuracy, 
within  the  current  state  of  the  art,  that  might  be  expected  when  predicting 
pressure  fluctuations  from  model  measurements.  One  obvious  factor  that 
influences  the  accuracy  of  flight  data  is  the  fact  that  flight-time  histories 
are  not  stationary. 

Since  the  RMS  level  of  amplitude  is  the  most  common  measurement  applied 
to  pressure  fluctuations,  it  is  appropriate  to  examine  the  effects  of  scaling 
on  these  measurements.  The  flight -time  histories  of  the  RMS  pressure  fluctua¬ 
tions  on  the  Mercury-Atlas  adapter  are  shown  in  figure  4.  Points  from  the 
7-percent  scale  model  extend  over  a  range  of  Mach  numbers  from  0.8  to  1.2. 
Since  there  were  differences  between  model  q0  and  full-scale  q0,  the  model 
results  have  been  adjusted  to  full  scale.  As  noted,  the  band-pass  frequency 
range  of  the  flight  data  extends  up  to  500  cps.  The  scaled  frequency  range  of 


*The  full-scale  data  were  obtained  from  unpublished  results  which  have 
been  compiled  by  Mr.  James  Ancell  of  Aerospace  Corporation. 
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the  model  should  extend  to  about  7,200  cps,  but  unfortunately  the  upper 
frequency  limit  of  the  tape  recording  of  data  was  only  6,000  cps.  Even 
though  the  band-pass  range  was  less  than  the  properly  scaled  range,  these 
7 -percent  model  data  still  generally  are  close  to  the  MA-5  data.  The  point 
indicated  by  an  x  shows  an  estimated  level  for  the  range  of  frequencies 
extrapolated  to  7,200  cps  from  the  power  spectrum  in  figure  3.  As  can  be 
seen,  the  agreement  between  the  7 -percent  model  data  and  the  average  flight 
data  is  improved.  Points  are  also  shown  in  figure  1  which  were  obtained  for 
a  band-pass  filter  range  from  3  to  5 00  cps  which  is  about  the  same  range  as 
for  the  full-scale  data.  The  difference  between  the  comparisons  of  model 
results  and  flight  data  for  the  two  filter  ranges  further  substantiates  the 
validity  of  frequency  scaling. 


Comparisons  of  Ranger  5  Results 


The  power  spectra  of  differential  pressure  fluctuations  on  the  Agena 
with  the  Ranger  5  payload  are  shown  in  figure  5*#  These  results,  which  were 
scaled  by  a  factor  of  about  10,  generally  agree  more  closely  than  did  the 
Mercury  data,  and  consequently  add  more  support  to  the  validity  of  scaling 
by  the  reduced  frequency  parameter.  It  should  be  noted  that  the  frequency 
range  of  the  measurements  of  fluctuations  on  the  Agena  was  very  low  compared 
to  the  range  available  on  the  Mercury.  This  fact  may  partly  account  for  the 
better  comparisons  between  the  Agena  measurements . 

Time  histories  of  the  RMS  levels  of  the  differential  pressure  fluctuations 
are  shown  in  figure  6.  As  with  the  Mercury  data,  the  range  of  band -pass  filter 
ing  of  the  signal  going  into  the  RMS  meter  was  scaled  inversely  as  the  model 
scale.  For  the  Agena  the  range  was  from  3  to  1,000  cps  for  the  model  and  from 
0.8  to  100  cps  for  full  scale.  With  these  band-pass  frequencies  precisely 
scaled,  the  maximum  intensities  were  reasonably  well  predicted  as  was  the 
variation  with  time  for  station  259-  At  station  219  the  flight -time  history 
indicates  an  earlier  buildup  of  intensity  and  also  a  dip  near  I7  seconds  that 
was  not  followed  by  the  wind-tunnel  data  which  were  taken  at  fixed  Mach  numbers 
One  might  speculate  that  these  differences  between  wind-tunnel  and  flight  data 
could  well  be  within  the  limits  of  accuracy  that  can  be  expected  considering 
the  unknown  effects  of  Reynolds  number  and  also  other  problems  associated  with 
obtaining  accurate  dynamic  measurements. 


*The  flight  power  spectra  and  a  magnetic  tape  of  the  differential  pressure 
fluctuations  were  obtained  from  Lockheed  Missiles  and  Space  Compamy.  The 
flight  data  will  appear  in  a  forthcoming  report  by  William  Eenricks,  Flight 
Test  Report  for  Ranger  Vehicle  6005,  SS/o26/5351>  I»'SSC/A384258. 
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Effects  of  Scale  on  Model  8 


The  separated  flow  on  the  Mercury  adapter  and  on  the  Agena  both  produced 
relatively  flat  power  spectral  densities.  Previous  experience  with  model  8 
for  other  investigations  has  shown  that  the  spectral  densities  ’/ary  in  shape 
from  a  predominately  low-frequency  spectrum  near  the  maxi  tram  diameter  shoulder 
to  a  flat  spectrum  as  the  distance  from  this  shoulder  is  increased.  For  this 
reason  model  8  was  originally  selected  for  the  investigation  of  scale  effects 
so  that  more  than  one  spectrum  shape  could  he  covered. 

Figure  7  shows  the  longitudinal  distribution  of  pressure  fluctuations 
on  model  8  for  three  different  model  sizes.  As  was  done  previously  the  band¬ 
pass  frequencies  for  the  RMS  readings  were  scaled  inversely  as  the  model  diam¬ 
eter  starting  from  an  upper  limit  of  2,000  cps  on  the  6.31-inch-diameter  model. 
Once  again  it  appears  that  satisfactory  agreement  was  obtained,  Further  sub¬ 
stantiating  the  use  of  the  reduced  frequency  for  scaling.  This  substantiation 
by  use  of  the  RMS  levels  is  only  valid,  however,  where  the  power  spectral 
levels  are  high  enough  in  the  frequency  range  being  scaled  to  effect  the  EMS. 
For  model  8  the  power  spectra  indicate  higher  energy  levels  at  the  higher  fre¬ 
quencies  at  stations  aft  of  x/'Dmx  =  1.4. 

As  previously  mentioned,  near  the  shoulder  on  model  8  the  pressure 
fluctuations  are  concentrated  at  the  lower  frequencies.  An  example  of  scaled 
power  spectra  from  a  station  within  this  region  is  shown  in  figure  8.  The 
results  appear  to  scale  well  including  the  peaks  at  wd/V  =  0.78,  which,  as 
a  result  of  the  satisfactory  scaling,  can  be  concluded  to  be  a  peak  caused  by 
an  aerodynamic  frequency  rather  than  a  model  motion.  It  is  also  interesting 
to  note  that  the  results  were  obtained  in  different  ’wind  tunnels.  The 
6.31-Inch  model  was  tested  in  the  6-  by  6-foot  tunnel  and  the  other  two 
models  were  tested  in  the  14-foot  wind  tunnel . 


Effects  of  Scale  on  Model  13a 


All  the  previous  data  shown  lead  to  the  conclusion  that  pressure - 
fluctuation  measurements  should  be  scaled  by  application  of  the  reduced 
frequency  parameter.  These  data  have  been  from  regions  of  separated  flow, 
and  the  spectra  have  been  smooth  curves  except  for  the  peaks  which  scaled 
on  model  8.  In  contrast  to  the  previous  smooth  spectra,  figure  9  shows  an 
example  of  power  spectra  of  fluctuations  in  the  region  of  the  shock  wave 
on  model  13a.  These  results  have  not  been  scaled.  It  can  be  noted  that 
several  peaks  coincide  on  the  two  smaller  models  which  were  sting  mounted, 
and  that  a  smoother  spectrum  without  such  predominant  peaks  was  measured  cn 
the  26-inch-diameter  half  model  which  was  mounted  on  the  tunnel  wall.  The 
fact  that  the  peaks  tend  to  coincide  when  the  model  supports  are  similar 
and  change  when  the  support  is  changed  indicates  that  the  shock-wave  motion 
is  influenced  by  model  motions.  The  peak  near  190  cps  also  coincides  with 
a  stream  disturbance  in  the  14-foot  wind  tunnel  (see  ref.  1)  thus  indicating 
an  influence  of  stream  fluctuations  as  well.  Since  it  would  be  expected  that 
the  full-scale  fluctuations  would  also  be  influenced  by  vehicle  motions  and 
stream  fluctuations,  the  details  of  spectra  of  this  type  certainly  could  not 
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be  scaled.  It  appears,  at  present,  that  for  design  purposes  the  best  approach 
would  be  to  construct  a  smooth  power  spectrum  having  a  slope  like  the  over -all 
slope  of  the  measured  curve  and  an  area  under  the  curve  equal  to  the  measured 
mean-square  amplitude.  There  is  some  justification  for  this  approach  since  as 
shown  in  reference -1  a  model  in  different  wind  tunnels  with  different  support 
properties  produced  the  same  RMS  measurements  even  though  in  one  case  a  large 
peak  predominated.  In  other  words,  it  appears  that  the  energy  available  from 
the  shock  wave  to  produce  fluctuations  is  unchanged,  but  that  the  frequency 
distribution  can  he  influenced. 


Effect  of  Total  Pressure  Variations 


The  power  spectra  of  the  pressure  fluctuations  have  been  put  in  coeffi¬ 
cient  form  by  dividing  the -spectra  measurements  by  q^-2.  While  generally 
this  has  been  assumed  to  be  a  correct  approach,  nevertheless,  some  measure¬ 
ments  were  made  on  model  8  to  check  the  effect  of  varying  q0  by  testing  at 
different  total  pressures.  Figure  10  shows  an  example  of  power-spectrum 
measurements  of  differential  pressure  fluctuations  obtained  in  the  11-  by 
11 -foot  transonic  wind  tunnel  at  total  pressures  of  from  15  to  60  inches  of 
mercury.  As  can  he  seen  there  was  little  effect  of  total  pressure. 


Correlation  of  Pressure  Fluctuations 


In  order  for  wind-tunnel  measurements  of  pressure  fluctuations  to  be  of 
value  it  is  necessary  that  the  spatial  correlation  of  the  fluctuations  on  the 
model  he  the  same  as  on  the  full-scale  vehicle.  There  has  been  insufficient 
time  to  obtain  cross  spectral  densities  of  any  of  the  scaling  data  before 
this  symposium;  however,  a  quick  look  at  the  over -all  correlation  was  obtained 
from  correlation  coefficient  measurements  as  shown  in  figure  11.  As  with  the 
RMS  values  previously  presented,  the  band-pass  frequency  range  of  the  signals 
to  the  analyzer  was  scaled  inversely  as  the  model  diameters.  While  there  are 
not  many  points  available  for  comparison,  those  shown  indicate  the  same  corre¬ 
lation  coefficients  for  a  difference  in  model  size  near  a  factor  of  1.  Other 
correlation  coefficient  measurements  with  respect  to  different  stations  show 
essentially  the  same  agreement  for  the  same  two  models. 


Attempts  were  also  made  to  perform  a  correlation  coefficient  analysis 
between  the  two  stations  tested  on  the  Ranger.  Average  coefficients  obtained 
from  the  flight  data  over  the  periods  of  time  from  1-3  to  46  seconds  after 
launch  and  from  1+5  to  50  seconds  agree  reasonably  well  with  the  10-percent 
scale  model  data  at  fixed  Mach  numbers  of  0.79  and  0.90,  respectively.  The 
corresponding  measurements  were  as  follows: 


Full  scale,  43-4-6  sec 
10-percent  model,  M  =  O.79 
Full  scale,  40-45  sec 
10 -percent  model,  M  =  0.90 


correlation  coefficient  =  0.155 
correlation  coefficient  =  0.130 
correlation  coefficient  =  0.100 
correlation  coefficient  =0.21 
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CONCLUDING  KEMAHKS 


On  the  basis  of  this  investigation  and  results  presented  herein,  it 
appears  that  pressure -fluctuation  measurements  on  models  in  regions  of 
separated  flov  should  be  scaled  to  full  scale  by  application  of  the  reduced 
frequency  parameter.  Scaling  by  a  factor  of  14. 3  on  the  Mercury  model  and 
by  a  factor  of  10  on  the  Ranger  5  model  gave  good  agreement  with  flight 
results. 

In  the  region  of  the  shock  wave,  model  measurements  can  be  influenced 
by  motions  of  the  model  and  also  by  disturbances  in  the  stream.  As  a  result, 
less  accuracy  can  be  expected.  It  will  probably  be  necessary  to  estimate 
power  spectral  densities  from  the  measured  RMS  level  and  a  generalized  spectrum 
shape. 

Correlation  coefficient  measurements  indicate  the  same  coefficients  for 
different  sized  models  of  model  8.  Reasonable  agreement  was  also  obtained 
between  the  10-percent  model  of  the  Ranger  5  and  full  scale;  however,  cress 
spectral  densities  are  required  before  definite  conclusions  can  be  drawn 
regarding  possible  effects  of  scale  on  the  spatial  correlation  of  pressure 
fluctuations . 
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Figure  1 . -  Models  investigated  for  scale  effects. 
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Figure  <?.-  Photograph  of  7 -percent  scale  Mercury  model. 
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Figure  3--  Power  spectra  of  pressure  fluctuations  on  the  Mercury -Atlas  adapter  near 


MA-4,  to  500  cps 
MA-5,  to  500  cps 
MA-7,  to  500  cps 
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Time  from  lift-off,  sec. 

Figure  4.-  Time  history  of  pressure  fluctuations  on  the  Mercury -Atlas  adapter. 
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Figure  5--  Power  spectra  of  differential  pressure  fluctuations  on  the  Agena  with  the  Ranger 

payload. 
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Figure  8.-  Power  spectra  of  pressure  fluctuations  on  model  8. 
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spectra  in  region  of  a  shock  wave 
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Figure  10 .  - 


Effect  of  total  pressure  on  power  spectra  of  differential 
pressure  fluctuations  on  model  8. 
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ABSTRACT 

Applications  of  flutter  models  in  the  design  of  high  performance  aircraft 
are  discussed  in  the  following  areas:  the  role  of  flutter  models,  scale 
factors,  design  and  fabrication,  support  system  requirements  and  stability  and 
control  for  complete  air  vehicle  models,  influence  coefficients  measurements, 
vibration  measurements,  flutter  test  instrumentation,  safety  devices,  tunnel 
flow  roughness  and  flutter  testing.  Information  is  presented  on  nev  construc¬ 
tion  materials  and  processes  Including  a  nickel-plating  technique  for  the 
application  of  external  skins,  &  novel  support  system  for  subsonic,  transonic 
and  supersonic  complete  air  vehicle  models,  a  novel  wind  tunnel  dynamic 
pressure  reduction  safety  device,  and  an  automatic  system  for  Influence 
coefficients  data  gathering  on  IBM  punched  cards. 


BOMEgCLATUEB 

b  «  reference  half-chord  length 
a>  =  circular  frequency 
V  m  true  velocity 

M  t=  mass  ratio,  M  ;  Poisson's  ratio,  G/e 
irpb2 

E  =  modulus  of  elasticity 
G  »  modulus  of  rigidity 

I  a  bending  moment  of  inertia;  mass  moment  of  inertia 
J  m  torsional  atiffhess  constant 
t  a  structural  thickness 
q  3  dynamic  pressure,  l/2  p  V2 

M  3  margin  built  into  a  flutter  model,  +/ 100;  Mach  number;  total  mass 
p  3  air  density 

ps  3  density  of  structural  material 
I  b  bending  moment  of  inertia  per  unit  length 
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Figure  Bo. 
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mass  par  unit  length 
static  deflection  or  deflection 
linear  spring  constant 
rotational  spring  constant 

area 

thickness  of  solid  aodel  i  it  nurture 
center  of  gravity 
elastic  center 
aerodynamic  center 


denotes  nodal 
denotes  full  scale 
refers  to  centnoidal  axis 
denotes  actual 
denotes  structure 
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IBTBOCTJCTIOI 

The  prevention  of  flutter  in  high  performance  aircraft  haa  become  a 
design  problem  of  continually  increasing  importance  as  flight  boundaries  have 
expanded.  This  has  been  so,  not  only  as  a  result  of  the  speed  increases,  but 
also  as  a  result  of  air  vehicle  configurational  changes  dictated  by  the 
broadened  operational  environment.  Inherent  stiffness  and  consequently  flutter 
speeds  have  tended  to  be  lower  as  structural  thickness  ratios  have  decreased, 
structural  temperature ■  have  Increased,  and  new  alloys  with  higher  strength  to 
stiffness  ratios  have  been  developed.  And  this  haa  not  been  completely  offset 
by  the  generally  favorable  effects  on  stlffbess  and  aerodynamic  parameters  of 
the  trend  toward  reduced  aspect  ratios  of  airfoil  surfaces.  The  net  result  of 
these  trends  has  been  a  tendency  for  Unit  speeds  to  approach  and  in  aurny  cases 
to  exceed  inherent  flutter  boundaries  resulting  from  strength  requirements  alor  -> 
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Prior  to  the  advent  of  supersonic  speeds,  designs  satisfying  strength 
requirements  vers  frequently  adequate  fro*  the  flutter  standpoint  except 
possibly  for  control  surface  aasa  balancing  and  approxlaate  flutter  analyses 
usually  sufficed  to  ensure  that  flutter  boundaries  eere  veil  In  excess  of 
Halt  speeds .  Vow  the  basic  structural  gages,  skin  panel  parameters, 
control  surface  stiffness  levels  and  Indeed  the  overall  configurations  are 
frequently  affected  by  flutter  considerations.  This,  together  with  the 
Increased  enphasls  on  design  optimisation,  the  Increased  cost  and  the  re¬ 
duced  susceptibility  of  sophisticated  high  performance  aircraft  structures 
to  quick  fixes  have  made  accurate  and  rapid  flutter  prediction  techniques 
an  economic  necessity.  Furthermore,  the  accurate  prediction  of  flutter 
boundaries  has  assumed  greater  importance  in  minimising  the  cost  and  risk 
of  flight  testing. 

However,  a  number  of  factors  have  combined  to  make  accurate  flutter 
predictions  more  difficult  and  tine  consuming:  (l).  The  trend  toward  low 
aspect-ratio  plate  like  structures  has  increased  the  difficulty  of  obtain¬ 
ing  accurate  analytical  determinations  of  structural  stiffness j  (2),  There 
has  been  a  tendency  for  design  evolutions  to  outstrip  the  development  of 
unsteady  aerodynamic  theories,  and  (3),  The  Increased  air  vehicle  complex¬ 
ity  and  the  greater  number  of  conditions  requiring  analysis  have  combined 
to  greatly  extend  the  time  and  effort  required  for  comprehensive  flutter 
evaluations. 

In  these  circumstances,  experimental  and  analytical  approaches  must 
be  carefully  combined  to  obtain  adequate  and  timely  solutions  far  flutter 
design  development  problems. 

The  present  paper  Is  limited  to  a  discussion  of  some  of  the  problems 
encountered  In  the  application  of  flutter  models  to  the  design  of  high 
performance  aircraft  and  techniques  ve  have  alloyed  in  their  solution. 


THE  BOLE  OF  FLPTT8R  lCBELS 


In  flutter,  as  In  other  fields  of  engineering,  experimental  and 
analytical  methods  complement  each  other  and  the  proper  balance  and  timing 
between  analytical  and  experimental  efforts  to  obtain  an  adequate  engineering 
answer  for  the  minimum  expenditure  must  be  determined  for  each  problem  on 
the  basis  of  prior  knowledge  and  experience  subject  to  revision  as  the 
problem  unfolds. 

Flutter  model  tests  have  the  following  specific  functions: 

1.  To  verify  analytical  natural  vibration  mode  and  flutter  predic¬ 
tion  methods. 

2.  To  evaluate  methods  far  computing  structural  stiffness. 
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3.  To  obtain  flutter  boundaries  by  direct  measurement 

4.  To  Investigate  parameter  trends 

Flutter  models  may  be  classified  In  tvo  broad  categories.  One  Includes 
the  simple  trend  models  and  the  other  the  detailed  models  of  an  actual  design. 
The  slaple  cantilever  trend  model  results  can  provide  approximate  or  prelimi¬ 
nary  Indications  of  Mach  number  and  other  parameter  trends  of  great  value  In 
formulating  a  more  comprehensive  attack  on  the  flutter  problem  and,  In  con¬ 
junction  with  detailed  complete  air  vehicle  analyses  end  low  speed  model  test 
results,  can  be  used  to  determine  flutter  boundaries  when  margins  are  not  too 
small.  When  greater  accuracy  is  required,  tests  of  detailed  component  models 
covering  the  range  of  Mach  numbers  can  be  used  to  supplant  the  trend  model 
results.  Of  course,  where  coupling  effects  between  components  are  not  small, 
adequate  margins  from  detailed  component  model  tests  become  reassuring  pre¬ 
liminary  results  which  must  be  supplemented  by  complete  air  vehicle  analyses, 
using  methods  which  correlate  with  model  results,  and/or  high  speed  tests  of 
complete  air  vehicle  models. 


SCA1E  FACTORS 


For  a  flutter  model  to  be  a  good  dynamic  representation  of  a  full  scale 
structure,  the  full  scale  stiffness  and  mass  distributions  and  the  external 
geometry  must  be  simulated  with  sufficient  accuracy.  In  addition,  there  should 
be  a  one-to-one  correspondence  of  bm,  M  and  M .  It  generally  Is  not  possible 

to  simulate  the  full  scale  Reynold^  number  but  the  effect  of  this  discrepancy 
on  flutter  probably  is  small,  except  for  Investigations  of  stall  flutter.  If 
model  Reynolds  numbers  are  above  the  critical  value. 

It  Is  Interesting  to  note  that  flutter  speeds  are  Independent  of  size. 

This  can  be  seen  from  the  facts  that  flutter  occurs  for  some  value  of  the  non- 
dimensional  parameter  bm  and  that  structural  frequencies  very  Inversely  with 

size  b,  thus  malxttalnl^  a  constant  value  of  bm  and  hence  flutter  velocity  ▼. 

There  are  basically  two  methods  for  obtaining  stiffness  simulation  In 
the  conversion  of  full  scale  structures  to  model  scale.  One  Is  the  direct 
scaling  of  the  full  scale  structure  and  the  other  Is  by  Indirect  scaling 
through  the  use  of  various  stiffness  coefficients.  The  method  of  direct 
sealing  Is  employed  for  structures  where  one -dimensional  stiffness  coeffi¬ 
cients  (EC,  GJ,  etc.)  are  not  applicable,  for  models  which  are  also  to  be  used 
to  define  or  check  full  scale  stiffness  and  for  high  speed  models  where  the 
correct  mass  ratios  must  be  simulated.  The  method  of  Indirect  scaling  Is 
employed  primarily  In  the  design  of  the  smaller  details  of  models  which  are 
-  otherwise  directly  scaled,  for  simple  research  models  and  far  lorn  speed  models 
which  can  be  described  In  terms  of  one -dimensional  stiffness  coefficients  and 
far  which  mass  ratio  simulation  can  be  achieved  with  structures  of  relatively 
low  stiffness  to  weight  ratios. 
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Is  general,  direct  scaling  tends  to  be  sore  accnrete  and  nodel  const ract ion 
ujr  be  accomplished  In  a  shorter  overall  tine  span  since  fever  analytical 
nanlpulations  of  the  structure  are  required. 


Sow  of  the  aost  useful  and  Important  xelatlonahlpa  between  flutter 
nodel  and  full  scale  structures  are  staeaarlsed  an  follows: 
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Model  structural  gages  may  be  computed  directly  from  full  scale  gages  using 
equations  (la)  and  Tib)  wherein  It  Is  assumed  that  structural  centroids  are 
matched  and  errors  In  Iq^/Iqj-  can  be  Minimized  sufficiently  by  proper  choice 
of  model  materials.  Where  -tee  full  scale  Poisson's  ratio  Is  not  exactly 
simulated  by  the  model  material,  better  accuracy  is  obtained  by  using  eq  (la) 
for  members  primarily  loaded  In  bending  and  eq  (lb)  for  members  primarily 
loaded  In  shear.  Relations  (2a)  and  (2b)  give  the  total  and  running  ratios  of 
model  to  full  scale  mass  respectively.  The  ratio  of  frequencies  is  specified 
by  relation  (3)  and  it  can  be  seen  that  the  reduced  frequency  parameter  is 
simulated  exactly  only  when  the  margin  M  built  into  the  model  is  zero. 

However,  it  is  common  practice  to  build  some  margin  (say  M  »  .15)  into  high 
speed  models  so  that  the  desired  stiffness  margin  can  be  demonstrated,  in 
tests  at  constant  Mach  number,  for  the  correct  mass  ratio.  The  ratio  of  static 
deflections  is  given  by  eq  (4).  Equations  (5),  (6),  (7)  and  (8)  give  the  ratio 
of  linear  spring  constant,  rotational  spring  constant,  torsion  constant  and 
moment  of  inertia  respectively.  Equation  (9)  specifies  the  area  ratio  for 
tension-compression  members  such  as  spar  caps,  longerons,  etc.  Equation  (10) 
can  be  used  to  calculate  the  thickness  of  solid  homogeneous  model  structures 
to  simulate  the  full  scale  El. 

It  is  interesting  to  note  that  fuel  slosh  frequencies  vary  inversely  as 
the  square  root  of  size  while  structural  frequencies  vary  Inversely  as  the 
first  power  of  size.  It  can  be  shown,  therefore,  that  slosh  frequencies  are 
properly  simulated  when  full  scale  static  deflections  are  properly  scaled. 


EESIQT  AHD  FABRICATION 


An  examination  of  the  fact  that  flutter  speeds  are  independent  of  size  can 
provide  some  insight  into  the  basic  problems  faced  in  the  design  and  fabrication 
of  flutter  models.  For  example,  if  we  were  to  conduct  model  teats  in  the  atmos¬ 
phere  using  free -flight  or  sled  techniques,  all  we  would  have  to  do  is  duplicate 
the  full  scale  structure  to  some  reduced  scale  in  all  important  details  using 
full  scale  materials  and  we  would  have  a  unity  speed  scale  flutter  model. 
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V. 


'  One  problem  that  night  he  encountered  is  that  the  aodel  structural  gages, 

1  having  been  scaled  according  to  the  geometrical  scale  factor,  would  he  too 

/  thin  to  fabricate.  We  might  solve  this  problem  by  using  a  material  vith  a 

lover  modulus  of  elasticity  requiring  correspondingly  greater  thickness  to 
provide  equivalent  stlfftaess.  But  the  substitute  material  mould  have  to  have 
a  stiffness  to  density  ratio  equal  to  or  greater  than  that  of  the  full-scale 
material  or  the  required  weight  could  not  be  net.  Also,  since  the  full-scale 
structures  are  designed  to  be  as  efficient  as  possible  In  terms  of  strength 
to  weight  ratio  and/or  stiffness  to  weight  ratio  and  the  materials  employed 
generally  have  a  high  stiffness  to  density  ratio,  lass  efficient  but  larger 
gage  concentrated  spar  types  of  model  construction  offer  no  possibility  of 
meeting  the  required  mass  ratio  unless  some  material  vith  an  ultra-high 
stiffness  to  density  ratio  such  as  beryllium  ia  used. 

Tor  models  designed  for  testing  in  wind  tunnels  the  problems  are  quite 
similar  and  may  be  stated  by  considering  the  pertinent  scale  factors.  The 
scale  factor  for  required  model  mass  is 

£  -  Aj if  £m  (2a) 

P  y 

and  the  scale  factor  for  actual  model  mass  is 

%  “  P*M  ^  (11) 

P  bf  *r  \h) 


Substituting  the  expression  for  tg/tj.  from  eq  (la)  into  eq  (ll)  gives 


Equating  actual  and  required  mass  scale  factors  gives 


Assuming,  for  exanple,  a  transonic  model  designed  for  a  blowdown  tunnel  using 
air,  the  ratio  of  Tu/7y  might  be  approximately  .95  and  with  a  built-in  margin 
of  15* 

(b/P,)m  *  (• 95 )2  1  =  .683  (1*0 
(s/P8)f  (1.15)2 
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This  npnmti  the  simian*  model  to  fall  actio  ratio  of  E/pg,  for  the 
assumed  conditions,  that  can  he  used  without  exceeding  the  repaired  weight. 
For  e  low  speed  Model  eq  (13)  with  M  -  0  can  he  solved  for  the  mart  mm 
velocity  scale  that  can  he  aohiered  for  given  structural  materials  without 
exeeedlag  the  repaired  weight  to  yield 


(15) 


Am  examination  of  eq  (la)  shows  that  to  maximise  model  gages  the  models 
should  he  as  large  as  possible,  he  built  of  materials  with  low  moduli  of 
elasticity,  he  tested  la  tunnels  capable  of  high  dynamic  pressures  and  not 
have  excessively  large  built -in  margins. 

In  summary,  model  construction  materials  with  low  elastic  moduli  and 
high  stlfAwsa  to  density  ratios  are  desired.  Table  I  lists  some  of  the 
Important  properties  of  materials  we  here  used  in  flutter  model  construction 
along  with  thosa  for  a  typical  full  scale  material  (steal)  for  oeaaparlion. 


TABLE  I 

MATERIAL  PROPERTIES 

9  9  pB  .  E/?8 


MATERIAL 

10°PSI 

10®PSI 

sU-  G/E 

LB/II.3 

lo^nr. 

Steel  (PH  15-7  Mo) 

29-5 

11.4 

•  39 

.28 

105 

Isowood  (23  lb/ft. 3) 

.0576 

.0221 

•  36 

.0133 

*.3* 

Isowood  (30  lb/ft.3) 

.0935 

•03*5 

•  37 

.0171 

5.48 

Styrofoam  (4.5  lb/ft.3) 

.00633 

.00225 

•  35 

.00211 

3.00 

Al.  Foil  (5052) 

10.3 

3.9 

.38 

.10 

103 

Vickel  (Electroplated) 

20* 

7.99* 

.*0 

.312 

64 

Al.  (7075) 

10.6 

4.0 

.38 

.10 

106 

Magnesium 

6.5 

2.* 

•  37 

.065 

100 

*  Higher  rallies  cam  he  obtained  by  varying  plating  bath. 
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Isowood  is  a  synthetic  isotropic  arterial  developed  and  Manufactured  to 
our  specifications  by  our  Space  sad  Information  Systems  Division  speci¬ 
fically  for  flutter  Model  applications.  This  Material,  due  to  its  low 
stiffness  and  low  stiffness  to  density  ratio,  is  suitable  primarily  for 
low  speed  Models.  for  high  speed  Models,  we  have  developed  a  technique 
of  electroplating  external  shins  of  nickel  over  a  Styrofoam  core  containing 
inlays  In  denser  Materials  representing  spars,  ribs,  bulkheads,  etc.  to 
obtain  highly  detailed  directly  scaled  Models.  The  advantages  of  the  nickel 
plating  technique  are  that  one  piece  skins  simulating  complex  tapers  and 
compound  curves  can  be  plated  and  to  lover  stiffness  or  St  values  than  are 
practical  with  the  aluminum  foil  technique. 

Some  examples  of  model  construction  are  shown  in  figures  1  through  6. 
The  low  speed  solid  Styrofoam  models  shown  in  figure  1  were  designed  to 
simulate  the  full  scale  SI  per  unit  length  throughout  the  planform  by  the 
application  of  eq  (10).  Many  small  lead  weights  are  imbedded  in  the 
Styrofoam  to  obtain  the  required  mass  distribution.  Because  of  the  short 
time  span  required  for  design  and  construction,  we  have  found  this  type 
of  model  to  be  very  suitable  for  preliminary  design  studies,  figure  2 
shows  some  detailed  skins  and  frames  machined  In  Isowood  and  Styrofoan 
for  directly  scaled  low  speed  flutter  Models.  In  some  cases  an  Internal 
sealant  is  applied  so  that  wet  fuel  can  be  used.  The  minimum  practical 
gage  in  Isowood  Is  approximately  0.020  Inches.  The  internal  structures  of 
some  directly  scaled  high  speed  models  are  shown  in  figure  3*  figure  4 
shows  a  directly  scaled  low  speed  complete  air  vehicle  model  fabricated 
largely  from  Isowood  and  Styrofoam.  Thin  sheet  as  well  as  block  types 
of  Styrofoam  were  vised  in  this  application.  High  speed  nickel  plated 
flutter  models  are  shown  in  figure  5»  In  figure  6,  the  uprper  portion 
shows  an  indirectly  scaled  low  speed  flutter  Model  and  the  lower  portion 
shows  some  of  the  external  wing-mounted  stores  constructed  for  It.  This 
model  was  constructed  of  many  separate  sections  with  each  supported  at  a 
single  point  on  a  concentrated  spar. 


SUPPORT  SI3IBC 

Rie  ideal  wind  tunnel  suspension  system  for  complete  air  vehicle 
flutter  models  may  be  defined  as  one  with  zero  mass,  stiffness  and 
damping,  capable  of  reacting  the  drag  and  normal  forces  without  creatix^ 
any  aerodynamic  disturbance  and  Including  very  accurate  automatic  model 
position  and  attitude  controls  which  will  not  affect  the  flutter  char¬ 
acteristics  to  be  determined.  Heedless  to  say,  this  is  not  the  system 
to  be  described  here.  However,  figure  7  illustrates  a  novel  suspension 
system  design  successfully  employed  in  low  speed  and  transonic  tests  of 
complete  air  vehicle  flutter  models  that  may  come  close  enough  to  the 
ideal  suspension  for  practical  purposes.  Utilizing  only  six  torque  bar 
springs,  the  system  provides  low  frequencies  in  all  rigid-body  modes 
except  fore-and-aft  translation  as  well  as  essentially  Independently 
adjustable  spring  constants  and  elastic  centers  for  vertical  translation, 
side  translation,  pitch,  yaw  and  roll. 
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In  addition;  the  effective  mu  and  duping  can  be  held  to  very  lew  values 
and  the  high  degree  of  aerodynamic  cleanliness  is  consistent  with  requires nts 
for  subsonic,  transonic  and  supersonic  testing.  When  required,  renote 
control  of  pitch  angle  can  be  Included  in  a.  straightforward  manner. 

Accurate  Indications  of  suspension  system  normal  force  and  pitching  moment 
reactions,  Indispensable  to  Intelligent  model  testing,  can  be  obtained  from 
strain  gage  bridges  applied  to  the  torque  bars  as  vlll  be  explained  later. 
Electrical  leads  from  model  strain  gages,  etc.  are  simply  routed  to  the  aft 
end  of  the  model,  bundled  together  and  run  through  a  sting  downstream  of  the 
model  with  a  suitable  amount  of  slack  between  sting  and  model.  Drag  loads 
are  carried  on  a  pair  of  cables  connected  to  a  whiffle  tree  upstream,  the 
whiffle  tree  is  supported  by  a  cable  sling.  A  pair  of  cables  with  coll 
sprites  la  used  to  preload  and  stabilize  the  vhif fie- tree -sling  assembly, 
these  springs  contribute  to  the  suspension  yaw  stiffness. 

Each  of  the  torque  bar  assemblies  is  identical  in  design  arrangement. 

A  typical  assembly  ia  illustrated  in  Figure  8.  the  outer  shell  is  a  stiff 
steel  tube  of  square  cross-section.  An  irreversible  worm  gear  drive  provides 
a  means  of  rotating  the  square  tube  about  its  longitudinal  center  line  in 
bearings  fixed  to  wind  tunnel  structure,  the  torque  bar  spring  is  a  solid 
steel  shaft  of  hexagonal  cross-section  mounted  concentrically  inside  of 
the  square  outer  shell  and  keyed  to  the  outer  shell  by  a  hollow  block 
machined  with  a  slip  fit  to  the  proper  shape,  the  keying  block  may  be 
moved  along  the  longitudinal  axis  of  the  spring  and  locked  in  place  by 
means  of  a  cable  and  pulley  system  and  a  series  of  set  screws  respectively 
to  change  the  torque  bar  effective  length.  A  sector  is  rigidly  attached 
to  an  adapter  collar  fitted  to  the  exposed  end  of  the  torque  bar  and  the 
adapter  collar  is  supported  in  ball  bearings  on  each  side  of  the  sector. 

A  cable  attachment  including  a  turbbuckle  is  fastened  to  the  sector  near  the 
hub  and  a  light  weight  cable  of  the  proper  length  is  routed  over  the  sector 
to  the  model.  A  swaged  ball  on  the  end  of  the  cable  seats  in  a  fitting 
built  into  the  model.  The  spring  constant  can  be  adjusted  by  sliding  the 
keying  block  to  change  the  torque  bar  effective  length  and  the  spring 
preload  can  be  changed  by  rotating  the  outer  shell  with  the  irreversible 
worm  gear. 

The  worm  gears  on  the  forward  and  aft  upper  and  lower  torque  bar 
assemblies  can  be  interconnected  by  shafting  and  appropriate  gearing  and 
driven  by  an  electric  motor  to  obtain  remotely  controlled  pitch  rotation 
about  an  axis  coincident  with  the  location  of  the  center  support  cables  as 
illustrated  in  Figure  7» 

The  cable  attach  fittings  in  the  model  can  be  designed  to  react  the 
rigging  preloads  in  direct  tension  and  hence  can  be  both  compact  and 
light  in  weight  and  their  weight  can  usually  be  completely  absorbed  as 
part  of  the  required  model  ballast. 
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For  aye trie  flutter,  suspension  system  effective  mass  it  ons -third 
of  tha  drag  cable  mass,  the  total  vertical  cable  aaei  aad  the  effective 
mate  of  the  aectora  aad  torqpe  bar  springs.  The  effective  seat  of  a  eeeter 
aad  torque  bar  tpriag  it 

**e  '  I  <* > 

12 


vheze  I  «  aatt  aotnt  of  iaertia  of  tector  pine  torque  bar 
tpriag  about  the  axle  of  rotation 

aad  r  «  tector  radiut 


Due  to  ltt  taall  radiut  of  gyration,  the  torque  bar  tpriag  aaket  a  negligible 
contribution  to  M*  aad  the  tector  contribution  caa  be  alaialKed  by  careful 
design.  For  tha  aatisynartric  cats  the  affective  aaat  of  the  suspension 
system  it  extremely  low  being  only  one-third  of  the  vertical  cable  neat, 
one -third  of  the  drag  cable  ease  for  lateral  translation  aad  the  total  drag 
cable  uses  plus  the  effective  anas  of  the  whiffle -tree  assembly  Including  the 
stabilising  cables  aad  springs,  in  year.  In  practice,  ratios  of  suspension 
effective  aass  (including  the  effective  mast  of  electrical  leads)  to  total 
model  mass  of  only  0.021,  0.003  sad  0.01^  have  been  obtained  for  pitch  and 
vertical  translation,  roll  aad  lateral  translation,  aad  yam  nodes  respectively. 

When  call  springs  are  used  to  suspend  a  model,  sene  oyronlse  on  spring 
stiffness  may  be  required  to  remove  the  coll  surge  frequencies  fna  the  flutter 
frequency  range  to  prevent  coll  surge  from  seriously  affecting  the  flutter 
characteristics.  This  problem  has  been  eliminated  through  the  use  of  torque 
bar  springs. 


Spring  constant  formulas  for  tha  symmetric  rigid-body  nodes  are  given 
below.  .  ^ 


f 
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Referring  to  figure  9 


\  "  \  +  \ 

%  “  %  +  \ 

%  “  *3b  +  *%,  ' 

the  spring  constant  In  re rt leal  translation  is 

St  -  Sv  4  %  ♦  Sv  * 

the  elastic  center  is  given  by 

£_  *27  ^12  4  K3Y  43 

^  Kly  +  Kgy  +  Kjy 


~\ 

> 


J 


(17) 


(IB) 


(19) 


and  the  pitch  stiffness  about  the  elastic  center  is  given  by 

°P  “  Kiv(sr)  +  Sv  ^  12  '  St)  +  *J9  ^13  "  St)  •  t20) 

Rearranging,  ee  obtain 


Iquations  (18),  (19),  end  (20)  can  be  solved  for  K-^y,  Kgy  and  K^y  ehen 
the  required  values  of  Kyj>,  Xyrj  and  Cp  are  given. 

Tor  the  purpose  of  illustrating  the  developawnt  of  spring  constant 
foraulas  for  the  aatisynuetric  rigid-body  nodes,  it  shall  be  assuned  for 
sinplidty  that  all  vertical  cable  lengths  are  equal. 
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Referring  to  figure  10,  the  following  set  of  equations  le  derived  for 
celculating  the  nodel  dead  weight  inactions  at  each  of  the  three  support 
points. 


wx  +  w2  +  w3  -  W  '  (a) 


The  differences  of  upper  and  lover  cable  tensions  are  given  fey 


where  the  upper  tensions  are  taken  positive  upward  end  the  lover  tensions 
are  taken  positive  downward. 


Referring  to  Figure  11,  the  lateral  stiffness  of  one  of  the  support  points 
Is  given  fey 


*L 

feut  6 

and 


Ta  sln9  +  Tl  sin  9 
8  6 
ie 

Tn  tl  *0  *  \ 

T  *  T  m  — 


for  small  deflections.  Then 


k1t  “  TlP  1  TlL 

L  i 

Kg^  »  T2U  +  T2L 


(25) 


and  the  total  spring  constant  in  side  translation  is 


=  kil  +  +  K3L 


(26) 
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The  vertical  location  of  tha  alaatlc  center  Is  given  by 


hf  1 1*  TJs  '  \  r0  *  l  (%  r%  '  \  \)  *  £  (\  %  ~  \  \) 


(27) 


% 

where  all  r's  axe  considered  positive  lengths. 
If  we  stipulate  that 


T  r 
U  ■  L 


(28) 


for  each  of  the  three  support  points,  each  of  the  bracketed  tens  in  eq.  (27)  is 
lero,  Zgp  >>  0  and  the  roll  axis  Is  coincident  with  the  X  axis.  Referring  to 
Figure  10,  the  fore-and-aft  location  of  the  elastic  center  is  given  by 

X  -  ^12  *  g3L  ^13  (29) 

“  *11  +  %  ♦  *3L 

and  the  yaw  stiffhess  about  the  elastic  center  is  given  by 

*  '  \  (vf  *  ■%  ({*  -  **)' 2  *  \  (*«  - 

Referring  to  Figure  12,  the  restoring  moment  in  roll  for  one  of  the  three 
supports,  assuming  small  deflections,  is 

Mr  *  ru6  U  +  TU  rUV  +  rL9  L  +  TjV  (31) 


(30) 


The  spring  constant  In  roll  is  given  by 
C 


.  *r  -  T_  r  4  .  9u)  »  t  r  4  * 
r  V  00  \  T/  ^  V  W 

Since  Gy  « ©L  «  y,  Kr  is  given  with  sufficient  accuracy  by 
Cr  "  TU  ru  +  ®L  rL 

Then  for  the  complete  suspension  system  the  roll  stiffness  is  given  by 
\  \  \)  *  (%  %  *  \  \)  *  (\  %  *  \  \ 


(32) 


(33) 
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Keferriag  to  Figure  13,  th«  contribution  of  the  whiffle  tro*  springs  to  y*» 

stiffness  la 

C ^  •  2  Ky  d2  (35) 


Stnnarlclng,  tha  aquations  for  tha  aatisyi—trlc  spring  constanta  ara 


*11  +  *2L  +  K3L  “  %r 


xst)  %,  +  '  ^12)  K2L  +  [*m  "  ^13)  k3L 


xst)  %  +ri2  ~  *sr)  *2L  +  ri3  "  k3l  "  ct  “  acvd 


cb  *  T: 


For  given  value  a  of  Kgj,  Xg^,  Cy  and  equations  (26),  (29)  and  (30)  can  be 
solved  for  the  required  values  of  K^,  and  K^.  Conblnlng  aq  (2*)  with 
eq  (25)  gives  the  following  expressions  for  the  cable  tensions: 
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Combining  relations  (28)  end  (3^)  gives 


CR  »  \  +  %  \  +  %  % 


(36) 


If  ve  assume  r,  =  r0  =  =  r  ,  then 


« 


CE  '  **  (%  *  %  *  % 


(37) 


or 


r”  ‘  ^ 7  %  *  *%) 


(38) 


and  using  relation  (28) 


r„  = 


r3L 


r  T10 
0  sr^ 
t1L 


rn  T2U 
u  «r— 

T2L 


r  ^32 

o  =r~ 

T3L 


*\ 


(39) 
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SUPPORT  SYSTPt  RgggrgBjBgS 


la  the  preceding  section  the  general  amngamt  of  a  novel  spring 
type  suspension  ayaten  la  described  la  detail.  However,  considerable 
engineering  is  required  to  neet  design  objectives  for  a  specific  application. 
These  requirements  Include  the  following: 

(1)  capability  of  adequately  supporting  the  dead  weight 

(2)  low  effective  mass 

(3)  low  damping 

(4)  low  stiffness 

( 3)  adequate  deflection  capability  to  handle  the  range  of  aodel  gross 
weights,  steady  air  load  reactions  and  dynamic  response  amplitudes 

(6)  sufficient  range  of  adjustment  on  spring  stlfftnss  and  cable  tension 
to  satisfy  static  and  dynamic  stability  requirements  for  all  loadings 

Some  flutter  models  are  so  fragile  that  they  can  fall  under  dead  weight  loads 
If  Inadequately  suppoi  -ed.  An  Important  feature  of  the  three  point  suspension 
system  is  that  it  can  be  adjusted  to  minimise  fuselage  vertical  bending 
deflections  and  stresses  due  to  dead  weight  loads.  Effective  mass  can  be  held 
to  a  minimum  by  designing  cables  and  sectors  to  low  strength  urgins  and  by 
specifying  minimum  requirements  for  deflection,  spring  constant  and  tension 
ranges.  Low  damping  is  assured  by  the  use  of  double  ball  bearings  to  support 
the  sectors  and  by  providing  adequately  large  cable  clearance  boles  In  the 
tunnel.  The  low  stiffness  requirement  tends  to  result  in  long  torque  bars. 

In  our  applications  the  torque  bars  were  18  to  20  feet  long.  The  required 
range  of  adjustment  of  spring  stiffness  and  cable  tensions  can  be  minimized  by 
optimizing  the  suspension  geometry  for  the  specific  application. 


SUPPORT  STSTOl  STABILITY  AID  00RTB0L 


Positive  static  or  weather  vane  stability  la  essential  to  the  maintenance 
of  low  trim  loads  and  Is  assured  as  long  as  the  support  system  elastic  centers 
In  pitch  and  yaw  are  forward  of  the  corresponding  aerodynamic  centers.  A 
finite  pitching  moment  coefficient  at  zero  lift  can  the  maintenance  of 
low  trim  loads  considerably  more  difficult.  The  pitching  moment  coefficient 
may  be  reduced  by  appropriate  control  surface  or  wing  deflection  If  the  local 
aerodynamic  loads  created  by  these  deflections  axe  not  excessive  and  the  effects 
on  aodel  stiffness  are  negligible.  Where  remote  control  of  aodel  Incidence 
angle  is  provided,  the  suspension  system  can  be  effectively  used  to  react  the 
pitching  moment. 
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The  qualitative  effect  of  center  of  gravity  on  rigid  body  longitudinal 
(or  lateral)  dynamic  stability  It  shorn  la  Figure  ik.  Regardless  of  the  pitch 
to  vertical  translation  frequency  ratio,  forward  centers  of  gravity  result  la 
rib  ratio*  sodas  vith  son-coalescent  frequencies  sad  hence  are  dynamically 
Stahl*.  Aft  C.G. 'a  are  show*  to  result  la  nodal  aerodynamic  stiffnesses 
teadlag  to  promote  frequency  coalescence  sad  dynamic  instehlllty. 

Our  experience  with  delta  winged  configurations  has  deaonstrmted  that 
dynsalc  Instehlllty  cat  be  obtained  for  a  forward  C.G.  with  >  1  a*  a 

«Vp 

result  of  the  anall  elastic  caster  deflections  In  the  vertical  translation 
node  that  can  occur  even  when  the  suspension  frequencies  are  more  than  two 
octaves  below  the  frequency  of  the  fuadaaental  elastic  node.  This  is 
Illustrated  In  Figure  15  for  a  two  point  support  system.  The  instability  was 
due  to  the  positive  Increment  In  aerodynamic  stiffness  In  the  vertical  trans¬ 
lation  nods  contributed  by  the  elastic  camber.  It  should  be  noted  that  camber 
deflections  are  shown  greatly  exaggerated  in  Figure  15  for  clarity.  For  a 
single  polst  support  It  is  shown  that  the  effect  of  camber  In  the  vertical 
translation  node  can  cause  divergence  when  °°p  >  1.  Fo r®p  >  1,  both 

®VT  ®VT 

problems  css  be  avoided  with  a  spring  constant  distribution  In  a  three  point 
suspension  representing  a  suitable  compromise  between  the  single  point  and 
two  point  extremes.  The  lower  part  of  Figure  15  Indicates  that  an  even 
better  solution  would  be  to  adjust  a  three  point  suspension  to  obtain  a  ratio 
of  <up  to  slightly  less  than  unity,  particularly  If  this  can  be  achieved 
In  a  specific  application  with  an  elastic  camber  deflection  contributing 
positive  aerodynamic  stiffhess. 


affmaras  oogffiuiKgTs  MBASPBBflnrrs 

The  measurement  of  model  flexibility  Influence  coefficients  Is  an 
extremely  important  aspect  of  flutter  model  testing.  These  data  are  used  to 
check  theoretical  Methods  of  stiffiaess  calculation,  to  assess  the  accuracy  of 
model  stiffness  against  calculated  or  neasured  full  scale  data  and  to  substan¬ 
tiate  the  accuracy  of  theoretical  methods  for  computing  vibration  nodes  and 
flutter  speeds  through  comparisons  with  measured  values. 

Unfortunately  these  measurements  can  be  both  tedious  and  tine  consualng. 

To  minimize  these  difficulties  and  at  the  saase  time  obtain  highly  detailed 
descriptions  of  model  flexibility  we  have  developed  the  use  of  automatic  data 
gathering  and  analysis  equipment.  The  heart  of  the  data  gathering  system  is  a 
hO  channel  Crescent  Engineering  nodal  5hQ  Commutated  Translator  unit.  Crescent 
variable  permeance  deflection  transducers  with  s  2  inch  stroke  providing  a 
maximum  overall  system  resolution  of  10.0002  inches  are  connected  to  the  inputs 
of  38  of  the  1*0  channels.  The  remaining  2  channels  receive  the  Inputs  from 
deflection  transducers  mounted  In  calibrated  force  rings  to  measure  the  applied 
loads.  The  force  rings  are  mounted  on  counterbalanced  arms  which  can  be 
rapidly  moved  and  locked  in  any  position.  Forces  are  then  applied  by  tumir.5 
coarse  screws  or  fine  differential  screws  Installed  between  the  force  rings  and 
the  supporting  arms.  The  electrical  outputs  of  the  force  rings  are  calibrated 
in  pounds  and  any  force  up  to  50  pounds  can  be  set  by  observing  the  output  on  a 
digital  voltmeter  in  the  Translator  unit. 
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The  da  flection  transducer*  consist  of  encased  center-tapped  colls  with 
central  probes.  Inch  probe  neighs  approximately  5  grans  sad  fits  loosely 
In  the  coll  assembly.  Typical  deflection  transducer  setups  consist  In 
■ousting  the  oolls  directly  below  the  nodsl  control  points  with  their  axes 
vertical  and  suspending  the  probes  on  strings  from  ns all  hooka  te^orarlly 
attached  at  the  nodal  control  points.  To  Measure  lateral  deflection#  or 
vhere  space  problems  exist  the  strings  are  routed  over  pulleys.  The 
Translator  unit  contains  an  automatic  printer  for  data  recording  and  an  IBM 
unit  has  also  been  connected  to  the  Translator  output  for  automatic  data 
recording  on  punched  cards.  Tbs  Translator  unit  contains  a  cross-bar 
scanner  ehich  automatically  sweeps  through  the  40  channels  In  approximately 
30  seconds.  In  this  Manner  individual  channel  sensitivities,  applied  loads 
and  deflections  are  automatically  punched  on  IBM  cards.  These  cards  are 
then  fed  into  a  computer  program  Milch  determines  the  best  slope  fit  of  the 
load-deflection  data  and  tabulates  these  slopes  a a  an  assembled  Influence 
coefficients  matrix.  Figure  16  shows  typical  setups  for  Influence  coeffi¬ 
cients  measurements.  With  this  equipment  the  data  for  a  30th  order  matrix 
have  been  recorded  In  less  than  two  hours  and  60th  order  matrices  are 
commonly  obtained. 


VIB  RATIO*  MEASUREMENTS 


The  measurement  of  natural  vibration  modes  is  also  a  very  important 
part  of  flutter  model  testing.  Our  first  natural  modes  measurements  on  a 
flutter  model  are  usually  made  in  the  Dynamics  Laboratory  for  the  boundary 
conditions  set  up  for  the  Influence  coefficients  measurements .  9>e  measured 

mode  shapes  and  frequencies  are  compared  with  calculated  values  baaed  an  the 
measured  influence  coefficients  to  obtain  a  check  on  the  accuracy  of  the 
measured  influence  coefficients  and  the  calculated  mass  data.  If  the  support 
conditions  In  the  wind  tunnel  are  different,  the  appropriate  analytical 
transformations  of  the  measured  Influence  coefficients  are  made  and  the  checks 
of  calculated  versus  measured  modes  are  repeated.  Figure  17  shows  a  typical 
multi -shaker  vibration  test  setup  In  the  wind  tunnel.  Good  agreement  is 
usually  obtained  between  measured  and  calculated  data  for  the  first  four  or 
five  modes.  Due  to  difficulties  In  accurately  measuring  higher  nodes  on 
lightweight  models,  more  reliance  is  usually  placed  an  the  calculated  higher 
nodes  and  the  calculated  modes  are  used  la  the  flutter  analyses  to  correlate 
theory  with  experiment.  In  complete  air  vehicle  flutter  analyses  we  have 
found  that  modes  as  high  as  the  ninth  free -free  mode  can  significantly  affect 
the  flutter  speed.  This  was  doe,  la  this  case,  to  an  overlap  in  the  frequency 
spectxums  of  predominantly  chordal se  fuselage  modes  and  predominantly  spamslse 
wing  modes. 
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flutter  aodel  wind  tunnel  instrumentation  should  be  designed  for 
simplicity  end  reliability  as  veil  as  for  el Blaue  weight  and  coat 
while  providing  the  following  Information  and  functions: 

(1)  Instantaneous  (pen)  recordings  of  aodel  response  waveforms 

(2)  Oscillograph  recordings  of  aodel  frequency  and  phase  response 
accurately  related  to  tunnel  test  conditions. 

(3)  Oscillograph  recordings  of  aodel  trie  and  dynamic  stresses. 

For  tests  in  continuous  tunnels  these  data  should  also  be 
displayed  on  aeters. 

(4)  Automatic  triggering  of  high  speed  action  picture  cameras  and 
aodel  safety  devices  for  tests  in  continuous  tunnels. 

In  our  experience,  strain  gage  sensors  have  proved  to  be  the  best  choice  for 
satisfying  the  above  requirements .  They  are  reliable,  light  in  weight  and 
low  In  cost.  Used  In  a  carrier  system,  their  ability  to  sense  frequency, 
phase  and  stress  and  the  possibility  of  feeding  more  than  one  circuit  with 
the  output  of  a  single  bridge  tends  to  promote  overall  simplicity. 

On  a  model  of  a  complete  air  vehicle,  six  or  more  strategically  located 
temperature  compensated  strain  gage  bridges  nagr  be  required.  For  the 
suspension  system  previously  described,  normal  force  Is  measured  by  installing 
a  complete  bridge  on  each  of  the  torque  bar  springs,  adjusting  the  outputs 
to  be  proportional  to  the  Individual  spring  constants  and  Interconnecting 
the  outputs  to  be  additive  for  pure  translation.  An  additional  set  of  six 
bridges  is  Installed  to  Measure  the  suspension  system  pitching  moment 
reaction.  A  total  output  proportional  to  pitching  moment  can  be  obtained 
by  adjusting  the  Individual  outputs  to  be  proportional  to  the  product  of 
spring  constant  and  distance  froa  the  elastic  center  and  staving  the 
outputs  with  proper  regard  to  signs. 


SAfBT  MYICBS 


A  quick  acting  tunnel  q  reduction  or  shutdown  device  is  required  If  one 
Is  to  plan  with  some  confidence  to  obtain  acre  than  one  flutter  point  from 
a  single  aodel. 

Figure  18  illustrates  one  of  the  earlier  devices  we  have  used  for  this 
purpose  In  our  low  speed  tunnel.  In  the  retracted  position  each  screen  Is 
held  flat  on  the  tunnel  floor  by  a  solenoid  operated  latch.  The  latch  Is 
actuated  by  an  electrical  signal  which  also  simultaneously  de-energises  the 
tunnel  fan  motor.  The  actuating  signal  can  be  applied  either  with  a  manual 
switch  or  automatically  froa  aodel  strain  gage  signals  through  an  electronic 
triggering  device  set  to  operate  at  predetermined  stress  levels.  When 
released  by  the  solenoid  operated  latch,  preloaded  springe  initiate  the 
upward  rotation  but  the  air  stress  supplies  most  of  the  force  for  extension. 
As  the  speed  brake  approaches  the  90  degree  travel  Halt,  hydraulic  shock 
struts  provide  a  smooth  deceleration  to  zero  velocity.  Full  extension  took 
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place  In  approximately  0.25  second*.  An  initial  q  reduction  of  vu 
obtained  over  the  area  blanketed  by  the  screen*  and  the  tunnel  Telocity 
decreased  to  zero  after  approximately  10  seconds,  The  screens  were  positioned 
in  the  tunnel  to  blanket  the  flow  orer  the  sings  of  complete  air  reticle 
■ode Is.  This  safety  derlce  was  first  used  with  the  model  shown  in  Figure  6. 
The  Initial  testa  sere  conducted  In  195**  and  approximately  1000  rune 
Including  400  flutter  points  hare  been  obtained  In  external  stores  flutter 
investigations  on  this  model.  This  model  Is  still  intact  largely  due  to 
the  effectiveness  of  this  "fly  swatter"  speed  brake  device. 

Figure  19  Illustrates  a  later  speed  brake  device.  It  consists  of  a 
false  floor  and  celling  with  trailing  edge  tabs  and  can  be  quickly  Installed 
In  the  test  section.  Like  the  earlier  device  It  is  coupled  to  the  fan 
motor  and  can  be  triggered  either  manually  froai  several  test  crew  stations 
or  automatically  at  predetermined  model  stress  levels.  Electric  motors 
release  the  tab  latches  and  drive  cams  to  initiate  tab  rotation.  The  air- 
stream  provides  most  of  the  energy  for  tab  rotation.  Rotation  through  the 
total  travel  of  100  degrees  occurs  In  less  than  0.1  seconds.  Figure  20 
shows  a  time  history  of  the  reduction  In  test  section  dynamic  pressure. 

The  additional  blockage  of  the  main  stream  and  the  aspirator  effect  on  the 
secondary  streams  due  to  tab  deflection  tend  to  divert  the  main  stream 
flow  to  the  secondary  streams.  In  addition  to  fast  response,  this  system 
provides  full  effectiveness  throughout  the  test  section  and  thus  can  be 
used  with  any  type  of  model  Installation. 


TWKL  FLOW  ROOGEHKESS 


While  some  tunnel  flow  roughness  or  turbulence  Is  desirable  for  the 
purpose  of  exciting  flutter,  it  can  create  serious  problems  In  the  testing 
of  flutter  models  with  low  mass  ratios  and  wing  loadings.  One  problem  Is 
that  the  apparent  flutter  speeds  obtained  may  be  somewhat  lower  than  the 
true  flutter  speeds.  A a  the  flutter  speed  is  approached  the  response  to 
turbulence  in  the  flutter  mode  increases  as  the  aerodynamic  daaplng 
decreases  with  the  result  that  the  true  flutter  speed  may  not  be  quite 
attained  due  to  (l),  misinterpretation  of  the  more-or-less  sustained  forced 
vibration  response  as  a  sustained  flutter  oscillation  and  (2),  the  desire 
to  avoid  overstressii^  the  model.  Another  possible  problem  Is  that  for 
even  lighter  and  more  fragile  models  the  stresses  Induced  by  flow  turbulence 
■ay  be  excessive  for  speeds  well  below  the  flutter  speed. 

Over  the  years,  no  flow  roughness  problems  were  encountered  In  tests 
of  relatively  high  mass  ratio  models  in  our  8x11  foot  low  speed  tunnel. 
However,  fairly  recent  attempts  to  test  low  mass  ratio  complete  air  vehicle 
flutter  models  indicated  that  the  turbulence  level  would  have  to  be  reduced 
if  these  tests  were  to  be  successfully  completed. 

Measurements  in  the  test  section  Indicated  small  fluctuations  in  both 
dynamic  pressure  and  flow  direction.  The  fluctuations  in  flow  direction 
were  considered  to  be  primarily  responsible  for  excessive  model  dynamic 
stresses  and  rigid  body  displacements— particularly  vfaen  the  testa  were 
conducted  at  zero  lift. 
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The  first  major  attempt  at  improving  the  flow  smoothness  vas  to  Install 
a  screen  vlth  about  50^  porosity.  This  vas  Installed  across  the  stilling  chamber 
just  upstream  of  the  contraction  to  minimise  power  consumption.  Subsequent 
flow  measurements  indicated  only  slight  Improvements .  After  nuch  additional 
thought  and  small  scale  experimentation  the  final  solution  vas  proposed. 

This  oonslsted  of  Installing  honeycomb  across  the  Mg  section  just  upstream 
of  the  screen.  The  boneyanfe  had  a  six  inch  square  cell  else  and  vas  three 
feet  long  in  the  streconrlse  direction.  Die  Improvement  vas  dramatic. 

Fluctuations  In  flow  direction  vere  reduced  to  between  1/5  and.  l/lO  of  the 
previous  level.  And  more  important.  In  teats  of  low  mass  ratio  flutter 
models,  the  onset  of  flutter  could  be  detected  on  a  pen  recorder  as  a 
relatively  clean  sinusoidal  response  often  before  any  motion  could  be 
detected  by  visual  observation  of  the  model. 

Later  on  when  the  same  problem  arose  In  attempts  to  teat  a  complete 
air  vehicle  model  at  transonic  speeds  In  our  7x7  foot  Trlaonlc  blow  down 
tunnel,  ve  knew  vhat  vas  required.  After  honeycomb  vas  installed  In  the 
stilling  chamber,  the  test  section  turbulence  vas  reduced  by  about  the 
seam  factor  as  for  the  low  speed  tunnel  and  the  transonic  complete  air 
vehicle  model  tests  vere  successfully  completed. 

PLOTTER  TBSTOG 

The  objectives  of  flutter  testing,  at  least  in  continuous  flow  tunnels 
and  particularly  for  lav  speed  models,  are  somevh m*  analogous  to  repeatedly 
valking  to  the  edge  of  an  abyss  and  safely  retret  I  ag.  It  Is  desired  to 
obtain  the  critical  flutter  speeds,  flutter  frequencies  and  mode  shape 
information  for  each  of  many  configurations  using  only  one  model.  This  is 
a  big  order,  and  Its  completion  requires  a  greet  deal  of  knowledge  and 
careful  planning  In  all  phases  of  the  woric  leading  up  to  and  including  the 
actual  testing. 

In  addition  to  the  preceding  portions  of  the  paper,  a  few  specific 
comments  may  be  In  order.  To  avoid  the  possibility  of  explosive  flutter, 
coulomb  friction  In  models  and  their  suspension  systems  should  be  held  to  an 
absolute  minimum,  sad  all  testing  should  be  conducted  veil  within  the  known 
model  strength.  An  Instantaneous  recorder  such  as  a  pen  type  is  extremely 
helpful  in  that  the  up-to-date  time  histories  of  flutter  burets  obtained 
make  possible  rapid  judgments  as  to  the  magnitude  of  damping  and  when  the 
point  of  sustained  flutter  has  been  reached. 

A  typical  pen  recording  of  a  flutter  point  obtained  on  a  low  speed  model 
is  shown  In  Figure  21.  Due  to  the  low  tunnel  flow  turbulence,  the  flutter 
response  is  extremely  clean  even  though  the  motion  is  so  small  as  to  be  barely 
perceptible  In  direct  viewing  of  the  model.  The  rapid  decay  of  the  flatter 
oscillations  following  the  actuation  of  the  "q  reducer"  also  may  be  noted. 

OOVCLODIKr  HBiAHCS 


It  is  clear  that  flutter  model  techniques  have  become  increasingly 
important  as  the  flutter  problem  has  assumed  greater  importance  in  the  design 
of  high  performance  aircraft.  In  the  future  ve  can  expect  further  developments 
in  dynamic  modeling  technology  as  new  problems  are  posed  by  more  advanced 
flight  vehicles. 
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Figure  1  Lov  speed  solid  styrofoam  flutter  models 
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Figure  4  Low  speed  directly  scaled  flutter  model 
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Figure  5  High  speed  directly  scaled  flutter  models 
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Figure  6  Lov  speed  indirectly  scaled  flutter  model 
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2  POINT  SUPPORT  fV/0  C,G‘  1  POINT  SUPPORT 
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Figure  15  Effects  of  model  elastic  camber  deflections  on  suspension  stability 
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High  speed  flutter  model 


Low  speed  flutter  model 


Figure  16  Influence  coefficients  measurements 
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Figure  18  An  early  wind  tunnel  speed  reduction  device 
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Tabs  extenaea 

Figure  19  Wind  tunnel  "q,  reducer"  safety  device 
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Figure  20  "q  reducer"  performance 
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Figure  21  A  typical  low  speed  flutter  record 
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ABSTRACT 


This  paper  discusses  the  design,  construction,  and  testing  of  a  dynamically  similar  model  used 
for  measuring  transient  loads  in  a  high  speed  wind  tunnel.  The  model  tested  during  this  program 
was  a  model  of  a  nuclear  ramjet  powered,  low  altitude  supersonic  vehicle  (LASV).  This  vehicle 
is  a  high  fineness  ratio  lifting  body  with  three  equally-spaced  stabilizers  and  a  canard  elevator 
system.  Dynamic  loads  at  two  stations  on  the  fuselage  obtained  for  store  ejection  and  transient 
elevator  inputs  are  compared  with  analytical  results. 
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L  INTRODUCTION 


The  content  of  this  paper  presents  data  obtained  under  contracts  AF33(657)-8686  and  AF33(657)- 
10439  with  the  United  States  Air  Force.  As  a  partial  fulfillment  of  these  contracts,  a  model  was  con¬ 
structed  and  tested  to  provide  empirical  data  for  comparison  with  analytical  studies.  The  objectives 
of  the  test  program  were: 

(1)  Assure  that  the  stabilizers  are  free  of  flutter  within  the  operational  limits  of  the  vehicle. 

(2)  Measure  transient  loadB  on  the  fuselage  resulting  from  the  ejection  of  an  internal  store. 

(3)  Measure  transient  loads  on  the  fuselage  resulting  from  elevator  inputs. 


IL  MODEL  DESIGN  AND  CONSTRUCTION 


Prior  to  the  finalization  of  the  model  design,  several  studies  of  the  model  were  required.  The 
size  of  the  model  and  the  values  of  various  nondimenslonal  parameters  had  to  be  established  and 
analyses  done  to  assure  the  suitability  of  the  values  chosen,  ft  was  decided  the  model  would  be 
built  without  the  inlet  and  diet  because  they  have  negligible  effect  on  the  symmetric  dynamical 
characteristics  of  the  vehicle. 

The  capabilities  of  the  vehicle  at  design  conditions  imposed  severe  requirements  for  dynamic 
simulation  with  a  model.  Because  of  its  size  and  dynamic  pressure  capabilities,  the  CVC  High 
Speed  Wind  Tunnel  was  chosen  for  this  series  of  tests.  This  tunnel  has  a  4'  x  4'  test  section  and  a 
dynamic  pressure  capability  at  Mach  3.0  of  about  4,800  psf.  The  large  test  section  permitted  the 
use  of  a  1/  20  scale  model  at  Mach  numbers  up  to  1. 6  without  interference  from  the  shock  wave  re¬ 
flected  from  the  tunnel  wall.  Also,  the  high  dynamic  pressure  capability  allowed  the  use  of  ?.  reason¬ 
able  stiffness  ratio. 
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The  model  geometric  scale  factor  is  defined  by: 

n-  1  -  0.05  U) 

m 


Where: 


n  is  the  geometric  scale  factor 

1^  is  a  representative  model  length 

1  is  the  corresponding  representative  vehicle  length 

a 

As  mentioned  above,  this  value  is  set  by  the  tunnel  size  and  lowest  Mach 
number  desired  for  testing. 

Hie  stiffness  ratio  of  the  model  is  determined  by  the  dynamic  pressure 
capabilities  of  the  wind  tunnel  used  for  testing.  The  dynamic  pressure  for 
the  vehicle  at  Mach  3-0  at  sea  level,  standard  MCA  conditions  is  13,346  pef, 
and  the  tunnel  is  capable  of  approximately  4,800  pef.  The  ratio  of  these  two 
dynamic  pressures  determines  the  stiffness  ratio.  It  was  decided  to  test  at 
a  design  dynamic  pressure  of  4-, 000  pef,  thus: 

f-^-0.30  (2) 

cja 


Where: 


f  is  the  stiffness  ratio 

is  the  dynamic  pressure  for  the  model 
q^  is  the  dynamic  pressure  for  the  vehicle 

The  mass  ratio  of  the  model  is  determined  by  the  air  density  capabilities 
of  the  wind  tunnel.  In  this  case  the  tunnel  is  capable  of  approximately  sea 
level  density,  therefore: 

■  -  Pm  =  1.0  (3) 

P  a 

Where: 

m  is  the  mass  raVo 

Pm  is  a  representative  model  density 

ps.  is  the  corresponding  vehicle  density 

The  model  frequency  scales  according  to  the  inverse  of  the  geometrical 
scale  factor  is  such  that  the  frequency  ratio  is  given  by: 
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r  »  ncu 

m 

Cl)  a 

Where: 

r  is  the  frequency  ratio 

is  a  representative  model  frequency 
^  a  is  the  corresponding  vehicle  frequency 


(4) 


The  frequency  ratio  is  related  to  the  stiffness  and  mass  ratios  of  equations 
(2)  and  (3)  respectively  such  that: 

r-j  ~f  ~  =  0.548 
J  m 


A  store  ejection  mechanism  was  designed  to  simulate  in  model  scale  the 
transient  loads  on  the  vehicle  resulting  from  the  ejection  of  an  internal 
store.  This  was  accomplished  with  a  pneumatic  device  that  imparts  a  scaled 
Impulse  to  a  brass  slug  of  scaled  mass. 

The  design  conditions  of  the  ejection  mechanism  were  to  eject  the  mass 
with  an  impulse  of  0.1148  lb^-sec.  The  force-time  history  of  the  ejection  was 
assumed  to  he  approximately  as  shown  in  the  non-dimensional  sketch  below. 


The  stroke  length  and  maximum  force  (F^...)  determined  from  this  curve  when 
0.0027  sec.  are  0.9  sec.  and  57. 40  lbf,  respectively. 

The  model  elevators  were  designed  to  introduce  transient  inputs  into 
the  system.  This  was  accomplished  with  a  slow  ramp  up  to  5°  deflection  of  the 
elevator  with  a  "quick  snap"  back  to  zero.  The  mass,  geometry,  and  root 
stiffness  of  the  elevators  were  scaled,  but  surface  flexibility  was  not. 

Studies  were  made  of  the  wind  tunnel  mounting  system  to  investigate 
the  effect  of  the  mounting  spring  stiffnesses  on  the  elastic  modes  and  fre¬ 
quencies  of  the  models.  The  coupled  system  was  theninvegtigated  to  insure 
that  the  system  would  be  stable  as  installed  in  the  wind  tunnel.  The  "free- 
free"  frequencies  and  the  coupled  frequencies  of  the  system  as  installed  on 
the  sting  are  shewn  in  Table  1.  The  elastic  mode  shapes  were  found  to  differ 
only  a  negligible  degree  from  the  free-free  to  the  constrained  condition. 

The  final  design  and  construction  of  the  models  were  subcontracted  to 
Dymanic  Devices,  Incorporated  of  Dayton,  Ohio.  The  primary  structure  of  the 
model  is  made  of  aluminum.  This  has  an  advantage  in  that  the  ratio  of  the 
Young's  modulus  of  aluminum  to  that  of  steel,  which  is  the  material  of  the 
vehicle,  is  almost  the  same  as  the  stiffness  ra+lo.  This  allows  many  pieces 
to  be  built  by  merely  scaling  the  vehicle  dimensions  down  to  model  size. 
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In  general,  the  philosophy  of  design  was  to  have  a  basic  structure  of 
proper  stiffness  with  a  light-weight,  non-load -carrying  aerodynamic  fairing. 

It  is  known  that  the  starting  loads  in  the  tunnel  can  be  as  high  as  16  pal 
across  the  model  at  M  =  3.0  with  7  or  8  pel  being  normal  at  this  Mach  nunfcer. 
These  loads  were  taken  into  consideration  during  the  model  design,  but  it 
was  not  possible  to  design  to  16  psi  and  still  maintain  the  scaled  stiffness 
required.  In  these  cases,  the  stiffness  design  was  used  in  the  hope  the 
resilience  of  the  model  mounting  would  alleviate  the  starting  loads  sufficien¬ 
tly  to  keep  from  damaging  the  model. 

A  general  arrangement  layout  of  the  fuselage  model  is  shown  in  Figure 
1.  The  structural  tube  was  attached  to  the  sting  by  means  of  leaf  springs 
designed  to  allow  pitch  and  translational  degrees  of  freedom  but  to  constrain 
the  model  in  yaw.  The  forward  mounting  was  designed  to  give  a  pinned  effect 
at  the  model  attachment  point  and  cantilevered  at  the  sting.  The  rear  support 
springs  are  clamped  at  both  ends,  both  at  the  model  and  at  the  sting.  The 
spring  attachment  points  were  located  as  near  to  the  node  points  of  the  first 
elastic  mode  as  possible  so  that  the  springs  would  have  a  small  effect  on  this 
mode.  Circular  steel  rings  were  used  as  stabilizer  mounting  reinforcing  rings. 
Brass  rings  were  attached  to  the  structural  tube  to  give  the  model  the  correct 
mass  distribution.  The  basic  structural  tube  and  attached  weights  were 
inclosed  with  a  balsa  shell  covered  with  silk  to  give  the  model  the  correct 
aerodynamic  contouring.  The  complete  model  and  a  cut-away  view  are  shown 
in  Figure  2. 

The  ejection  mechanism  is  a  pneumatically-operated  piston-cylinder 
arrangement.  The  slug  (which  is  the  piston  in  this  case)  is  held  in  place 
by  a  brass  screw  machined  to  the  proper  diameter  to  break  at  57  pounds  load. 
The  calculations  indicate  that  this  breaking  load  should  give  an  impulse  of 
about  0.11  lbf-sec.  Air  pressure  supplied  to  the  cylinder  provided  the  force 
required  to  break  this  screw.  The  ejection  mechanism  was  positioned  on  the 
antinode  of  the  first  elastic  mode  to  insure  the  maximum  response  in  this 
mode.  The  location  of  the  ejection  mechanism  is  shown  on  the  general  ar¬ 
rangement  drawing  of  Figure  1. 

Movable  elevators  made  of  solid  aluminum  were  placed  on  the  model  so 
that  transient  inputs  could  be  induced  upon  the  model  during  the  tunnel  runs . 
The  elevators  were  driven  by  an  electric  motor  through  a  cam  arrangement 
which  lifted  the  elevators  to  an  angle  of  5°  and  allowed  them  to  snap  back  to 
zero  through  a  half  cycle  of  the  cam.  The  elevators  were  returned  to  the 
zero  position  with  a  spring  arrangement.  The  electric  motor  was  mounted  at 
the  fuselage  bulkhead  at  the  cone-cylinder  break  (this  is  also  the  point  of 
attachment  of  the  forward  support  springs ) .  The  location  of  the  elevators 
and  the  actuator  motor  is  shown  on  Figure  1. 

The  initial  model  stabilizer  construction  was  basically  the  same  as 
the  vehicle  stabilizer  structure  used  as  the  basis  of  the  mathematical 
model.  A  planform  drawing  of  the  stabilizer  model  is  shown  in  Figure  3.  The 
stabilizer  structure  was  later  modified  to  the  planform  shown  in  Figure  h  to 
give  the  model  added  strength.  The  initial  stabilizer  model  was  made  up  of 
three  beams  with  the  space  between  the  beams  filled  with  balsa.  The  beams  and 
balsa  were  covered  with  an  aluminum  skin  which  was  brought  up  to  proper  aero¬ 
dynamic  contour  with  balsa  and  covered  with  Japanese  tissue.  The  distribution 
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of  corrective  masses  for  the  model  stabilizers  was  determined  from  the  data 
obtained  by  cutting  up  a  pilot  moCul  without  corrective  ballast  and  measuring 
its  mass  distribution.  The  difference  between  this  distribution  was  made  up 
with  thin  lead  sheets  attached  to  the  aluminum  skin.  The  later  design  was 
made  by  machining  the  beam  grid  from  a  solid  plate  and  filling  the  voids  with 
balsa.  The  corrective  masses  were  then  added  and  the  stabilizer  was  cover e' 
with  Japanese  tissue.  IMs  construction  was  found  to  yield  a  satisfactory 
simulation  of  the  plate-like  vehicle  stabilizer. 

A  sting  mount  was  designed  for  the  purpose  of  running  flutter  tests  with 
a  single  stabilizer.  Root  attachment  springs  were  included  to  simulate  the 
root  stiffness  conditions  of  the  fuselage. 

All  components  of  the  model  were  thoroughly  checked  by  CVC  personnel 
prior  to  acceptance  of  the  model  from  the  model  subcontractor.  The  stiffness 
and  mass  properties  of  the  models  were  closely  checked  and  modified  where 
necessary  in  an  attempt  to  keep  these  properties  to  within  15$  of  the  design 
values.  However,  some  compromises  with  stiffness  had  to  be  made  in  the  interest 
of  greater  strength. 

The  stiffness  of  the  fuselage  structural  tube  was  measured  in  the 
following  manner.  The  tube  was  cantilevered  at  the  tail  cone  break  point  and 
loaded  with  a  shear  load  at  the  location  of  the  elevators.  Mirrors  were  placed 
along  the  tube  at  equal  intervals  and  the  angular  deflection  of  light  beams 
reflected  from  these  mirrors  measured  for  the  loaded  condition.  These  deflec¬ 
tions  were  reduced  to  stiffness  distributions  and  compared  with  scaled  vehicle 
stiffness  curves  and  calculated  model  stiffness  curves.  This  comparison  is 
shown  in  Figure  5.  A  plot  of  the  distribution  of  the  stiffness  ratio  along 
the  fuselage  is  shown  on  the  same  figure  based  on  the  calculated  stiffness 
distribution.  The  model  stiffness  was  high  in  the  midsection  because  of  the 
difficulty  in  machining  the  fuselage  tube  to  within  less  than  0.010-inch  wall 
thickness . 

The  inertial  properties  of  the  completely  assembled  model  were  measured 
as  follows.  First,  the  model  mass  wa s  measured  with  an  accurate  balance  and 
the  center  of  mass  located.  The  moment  of  inertia  about  the  center  of  gravity 
was  then  measured  by  swinging  the  model  as  a  bifilar  pendulum  with  the  axis 
of  rotation  through  the  center  of  mass.  The  comparison  of  model  and  design 
mass  data  are  shown  in  Table  2. 

A  pendulum  type  device  was  made  to  measure  the  impulse  of  the  ejection 
mechanism.  Hie  simulated  store  was  fired  horizontally  from  the  ejector  mounted 
on  the  pendulum  and  the  maximum  amplitude  of  the  pendulum  swing  measured .  This 
amplitude  was  related  to  the  impulse  of  the  ejection  mechanism  by  the  use  of  the 
geometry  and  inertial  characteristics  of  the  pendulum.  After  the  ejection 
mechanism  was  calibrated  to  yield  approximately  0.11  lbf-sec  impulse  consistent¬ 
ly,  it  was  installed  in  the  fuselage  tube.  The  ejection  mechanism  was  static 
fired  several  times  prior  to  the  tunnel  runs,  and  high-speed  photographs  were 
taken  of  an  ejection.  The  impulse  was  calculated  by  measuring  the  exit  velocity 
from  these  sequence  photographs  and  multiplying  it  by  the  mass  of  the  slug. 

This  resulted  in  a  precise  measurement  of  the  ejection  impulse  of  0. 1095  lbf- 
sec.  The  structural  damping  was  calculated  from  a  trace  showing  the  bending 
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moment  time  history  at  a  point  18.38  inches  from  the  nose.  The  structural 
damping  was  found  to  he  2 %  of  the  equivalent  critical  viscous  (linear) 
damping. 


Initial  testing  of  the  stabilizers  was  done  with  a  mounting  block 
designed  to  simulate  a  cantilever  condition  at  the  stabilizer  attachment  points. 
This  block  was  clamped  to  a  rigid  test  fixture  and  mirrors  placed  on  the 
stabilizer  at  the  panel  points.  The  stabilizer  was  loaded  and  angular  deflec¬ 
tions  of  these  mirrors  in  the  x  and  y  directions  measured  by  the  light  beam 
method  mentioned  above.  These  measured  angular  deflections  were  compared  with 
corresponding  data  derived  analytically.  A  sample  of  these  comparisons  is 
shown  in  Figures  6  and  7  for  the  x  and  y  angular  deflections  respectively. 

The  mass  of  each  completed  stabilizer  was  measured  with  an  accurate 
balance  and  the  center  of  mass  was  located  by  balancing  it  along  two  axes  with 
a  knife  edge.  The  moments  of  inertia  of  each  stabilizer  were  measured  about 
two  mutually  perpendicular  axes  by  swinging  each  stabilizer  as  a  bifilar 
pendulum.  As  mentioned  earlier,  a  pilot  model  was  constructed  and  cut  up  for 
the  purpose  of  obtaining  mass  distributions  of  the  models.  Each  individual 
bay  was  weighed  and  the  corrective  masses  and  locations  thereof  required  to 
bring  the  models  into  agreement  with  the  design  values  were  calculated.  Lead 
weights  were  added  to  each  subsequent  model  to  provide  these  corrections. 

The  root  stiffness  of  the  sting  root  mounting  springs  was  measured  with 
a  combination  of  a  mirror  to  detect  angular  motion  and  a  dial  gage  to  detect 
lineal  motion.  The  comparison  between  the  design  root  influence  coefficient 
matrix  and  the  measured  root  influence  coefficient  matrices  are  shown  in 
Figure  8. 

An  acoustic  shaker  was  used  to  excite  the  stabilizer  models  for  the 
natural  frequency  checks.  A  vibration  pickup  with  its  output  connected  to  an 
oscilloscope  was  used  to  detect  the  motion  of  the  models  and  salt  was  sprinkled 
over  the  surface  to  define  the  node  lines.  A  sample  comparison  of  the  measured 
and  calculated  node  lines  and  frequencies  for  the  cantilevered  root  condition 
are  shown  in  Figure  9. 


III.  TESTING  TECHNIQUES 


All  wind  tunnel  tests  were  conducted  in  the  CVC  High  Speed  Wind  Tunnel. 
The  characteristics  of  this  tunnel  are  summarized  in  Figure  10.  This  tunnel 
was  chosen  because  of  its  large  size  ( h- '  by  4'),  wide  Mach  number  capability 
(0.2  to  5.0),  and  high  dynamic  pressure  capability  (5,000  pef  at  Mach  3*6). 

The  disadvantage  of  an  unevacuated  blowdown  tunnel  is  the  existence  of 
high  starting  loads  when  running  supersonic  tests.  These  loads  resulted  in 
several  premature  structural  failures  of  both  the  stabilizer  and  the  fuselage. 

A  stabilizer  was  mounted  on  a  sting  with  root  supports  which  simulated 
the  bulkhead  flexibility  for  the  purpose  of  flutter  testing.  The  frequencies 
of  each  model  being  tested  were  checked  after  each  run  to  check  for  possible 
structural  failures  during  the  run.  This  served  the  dual  purpose  of  checking 
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for  flutter  of  the  surface  and  indicating  the  serverity  of  the  tunnel  starting 
and  unstarting  loads. 

Each  of  the  stabilizers  was  instrumented  with  two  sets  of  strain  gages. 

One  set  was  oriented  to  pickup  s panvise  loads,  and  the  other  set  was  oriented 
to  pick  up  chordvlse  loads.  Each  of  the  strain  gage  traces  was  recorded  on  an 
oscillograph  type  recorder  simultaneously  with  several  tunnel  parameters.  In 
addition,  each  test  was  covered  with  a  high-speed  movie  camera  running  at 
approximately  1,000  frames  a  second.  A  stabilizer  was  successfully  flutter 
tested  up  to  Mach  2.0.  The  model  was  stable  during  the  run  but  was  broken  by 
the  unstart  shock. 

A  set  of  stabilizers,  modified  for  greater  strength,  was  fabricated  for  a 
second  series  of  tests.  The  most  important  modification  was  that  the  strength 
of  the  rear  beam  was  doubled.  It  was  found  that  the  beam  was  then  strong 
enough,  but  that  the  skin  (which  is  only  0.002-inch  thick  2024 -ST  aluminum) 
began  breaking  directly  behind  the  middle  beam  on  high  Mach  number  runs.  Thus, 
although  one  weak  point  was  eliminated  another  was  brought  to  light. 

The  complete  dynamic  model  was  mounted  on  the  sting  for  tunnel  testing  as 
shown  in  Figure  2.  The  balsa  fairing  was  broken  on  the  first  run  of  the  first 
series  of  tests  necessitating  that  the  model  be  modified  and  rebuilt  before 
testing  could  continue.  Even  after  the  model  had  been  strengthened,  the  only 
data  obtained  without  using  a  protective  device  was  at  Mach  1.6. 

The  response  due  to  elevator  deflections  was  measured  first.  The  elevator 
actuator  motor  was  controlled  by  a  switch  in  the  tunnel  control  room.  The 
elevator  was  deflected  slowly  up  to  5°  and  snapped  back  with  a  spring.  After 
the  shock  had  gone  down  the  tunnel,  the  elevators  were  activated  and  several 
responses  for  elevator  inputs  were  measured.  Next,  the  response  due  to  a  store 
ejection  was  measured.  The  air  pressure  value  was  controlled  by  a  switch  from 
within  the  tunnel  control  room.  The  shape  of  the  force-time  curve  has  been 
described  previously. 

The  model  was  badly  damaged  with  the  starting  and  unstarting  shocks  on  the 
first  run  at  Mach  2.2.  The  left  stabilizer  was  lost  during  the  start,  but  the 
right  stabilizer  remained  intact  during  the  run  and  apparently  was  lost  during 
the  tunnel  shutdown.  This  ended  the  attempts  to  test  this  model  without  some 
sort  of  protective  device  to  shield  the  model  from  the  high  loads  experienced 
during  tunnel  start  and  unstart. 

A  set  of  proximity  plates  were  designed  and  installed  in  the  tunnel  to 
protect  the  model  during  the  starting  and  unstarting  conditions.  The  plates 
consisted  of  steel  channels  which  completely  enclosed  the  model  on  all  sides 
and  were  open  on  both  ends.  The  plates  were  hydraulically  operated  on  command 
from  the  tunnel  control  room,  requiring  approximately  5  seconds  for  actuation. 
The  model  was  enclosed  during  the  tunnel  start  and  then  the  plates  were  retrac¬ 
ted  to  the  tunnel  wall  while  data  were  measured.  When  the  desired  information 
had  been  recorded,  the  plates  were  extended  to  protect  the  model  again  during 
the  unstart  conditions.  It  is  estimated  that  the  protective  plates  reduce  the 
starting  and  unstarting  loads  approximately  50$  to  75$  in  the  Mach  2. 0-3.0 
range.  However,  high-oscillating  loads  were  introduced  while  the  model  was 
fully  enclosed  by  the  plates,  with  Immediate  relief  of  these  loads  noted  upon 
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initial  plate  movement.  Separating  the  plates  initially,  approximately  l/U 

inch  from  the  closed  position,  resulted  in  considerable  reduction  of  these 

oscillating  loads.  This  method  of  protecting  models  is  in  a  developmental 

stage  and  requires  further  study  and  tests.  * 

Strain  gage  instrumentation  vaa  used  at  tvo  stations  on  the  model,  at 
21.3  inches  and  28.8  inches  from  the  nose.  The  strain  gages  vere  mounted  * 

directly  on  the  top  and  bottom  of  the  fuselage  structural  tube.  The  instru¬ 
mentation  was  calibrated  by  noting  the  deflections  of  the  oscillograph  traces 
for  known  loading  conditions  and  obtaining  a  scale  factor  in  inch-pounds  per 
inch  of  trace  deflection. 

Immediately  prior  to  each  run,  the  model  was  shaken  with  an  electromagnetic 
shaker  to  determine  the  first  few  natural  frequencies.  This  was  done  as  a  check 
on  the  structural  soundness  of  the  model  on  the  premise  that  any  failures  would 
manifest  themselves  in  frequency  shifts  in  the  model. 

Finally,  the  proximity  plates  were  closed  to  protect  the  model  during  the 
tunnel  start  and  the  actual  wind  tunnel  test  begun.  After  the  starting  shock 
had  gone  down  the  tunnel,  the  proximity  plates  were  retracted  and  the  high-speed 
movie  cameras  started.  The  store  ejection  was  automatically  triggered  to  occur 
1  second  after  the  start  of  the  movie  cameras  and  the  elevators  were  actuated 
manually  approximately  1  second  after  the  store  ejection.  After  a  few  readings 
had  been  obtained  for  elevator  inputs,  the  proximity  plates  were  replaced  to 
protect  the  model  from  the  unstart  shock  and  the  tunnel  was  shut  down.  Each 
wind  tunnel  run  required  approximately  15  seconds  of  run-time. 

The  data  points  obtained  in  this  manner  are  summarized  in  Table  3.  ' 


IV.  TEST  RESULTS 


Vertical  bending  moment-time  histories  were  obtained  for  both  elevator  and 
6tore  ejection  inputs  at  the  conditions  summarized  in  Table  3. 

It  was  found  that  the  model  was  bottoming  out  against  the  sting  for  the 
initial  elevator  deflection,  thereby  complicating  the  interpretation  of  the 
results.  However,  the  data  obtained  was  quite  clear  with  a  minimum  of  back¬ 
ground  "noise"  from  the  tunnel  turbulence.  A  reproduction  of  a  typical  data 
trace  for  elevator  inputs  is  shown  in  Figure  11. 

The  data  obtained  for  the  store  ejections  had  no  extraneous  inputs  as 
encountered  in  the  elevator  deflection  tests,  lhe  data  obtained  wus  very  clear 
for  all  conditions.  A  reproduction  of  a  typical  data  trace  for  store  ejection 
is  shewn  in  Figure  12  and  a  time  sequence  of  the  store  trajectory  at  M  =  3.0  is 
shown  in  Figure  13 . 

The  model  construction  and  testing  are  discussed  extensively  in  reference  1. 
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V.  CCMPARISCHS  WITH  ANALYTICAL  STUDIES 


<  In  defining  the  physical  system  and  formulating  the  mathematical  model  for 
use  with  a  digital  computer,  the  continuous  system  is  replaced  with  a  discrete 
system  of  collo&ation  points  by  the  application  of  certain  interpolation  proce- 

<  duxes  and  Lagrange's  equations  of  motion.  The  model  fuselage  was  idealized  as 
a  one-dimensional  body  whose  aeroelastic  properties  are  described  as  functions 
of  a  single  special  coordinate  along  the  fuselage.  The  first  four  calculated 
mode  shapes  and  frequencies  are  shown  in  Figure  14  compared  with  measured  data 
from  the  model  test. 

For  the  purposes  of  analysis,  the  stabilizers  are  idealized  as  plates  with 
integral  stiffeners.  The  aeroelastic  properties  are  described  as  functions  of 
two  spacial  coordinates  in  the  plane  of  the  stabilizer.  The  first  six  calculated 
node  lines  and  frequencies  cure  compared  with  measured  node  lines  and  frequencies 
for  a  typical  stabilizer  with  flexible  supports  in  Figure  15. 

The  results  of  the  "bent  body"  pressure  distribution  measurements  of  ref¬ 
erence  2  were  used  as  the  basis  of  the  development  of  the  aerodynamic  influence 
coefficient  matrices  for  Mach  3*0  flight  conditions.  Two-dimensional  piston 
theory  was  used  for  the  aerodynamics  of  the  stabilizers  and  elevators. 

The  methods  of  analysis  are  discussed  in  detail  in  references  1  and  3. 


VI.  COMPARISONS  OF  TEST  RESULTS  AND  ANALYSIS 

* 


The  primary  purpose  of  these  tests  was  to  provide  an  empirical  check  on  the 
assumptions  and  methods  utilized  to  calculate  the  transient  loads  on  the  vehicle. 
Because  it  was  not  feasible  to  achieve  a  free -flight  condition  in  the  wind 
tunnel,  the  model,  including  the  tunnel  mounting  system,  was  analyzed  for  direct 
comparisons.  The  philosophy  here  is  to  approach  as  complete  a  simulation  as 
possible  with  the  model  and  assume  that  the  assumptions  and  method  of  analysis 
are  adequate  for  the  full-scale  vehicle  if  they  are  adequate  for  the  model. 

Transient  responses  of  the  mathematical  model  were  calculated  by  solving 
the  equations  of  motion  with  a  numerical  integration  technique.  The  force  in¬ 
put  time-history  used  for  the  store  ejection  analysis  is  the  same  as  the  one 
described  in  Section  II.  The  elevator  input  analysis  was  approached  as  an  ini¬ 
tial  deflection  problem.  It  was  assumed  that  the  ramp  up  was  slow  enough  to 
induce  insignificant  dynamic  response  in  the  model  until  the  elevators  were 
snapped  back  to  zero  with  the  return  spring  (see  Figure  11 ) .  Thus,  the  initial 
deformed  shape  for  a  static  elevator  deflection  of  5°  was  calculated  and  used 
as  the  static  initial  condition  for  an  analysis  of  the  dynamic  system  with  zero 
elevator  deflection.  As  mentioned  earlier,  it  was  found  that  the  model  was 
bottoming  against  the  sting  initially,  requiring  that  the  problem  be  run  in 
two  parts  to  obtain  the  complete  solution.  Dynamic  loads  were  calculated  from 
these  responses  by  application  of  the  modal  acceleration  concept  for  calcula¬ 
ting  collocation  point  loads. 


145 


CONFIDENTIAL 


CONFIDENTIAL 

A  comparison  of  the  calculated  and  measured  bending  munents  for  a  zero 
airspeed  store  ejection  is  shewn  in  Figure  16  and  comparisons  for  store  ejec¬ 
tion  and  elevator  Inputs  at  Mach  3  conditions  are  shown  in  Figures  17  and  18 
respectively.  As  a  matter  of  Interest,  the  calculated  load-time  history  for 
a  store  ejection  on  a  system  without  the  sting  support  springs  (free-free)  is 
included  in  Figure  17. 


VII.  CONCLUSIONS 


It  has  been  demonstrated  that  relatively  fragile  dynamic  models  can  be 
tested  successfully  at  supersonic  speeds  in  unevacuated  blow-down  wind  tunnels 
if  protective  devices  are  employed  to  reduce  the  start  and  unstart  loads  an 
the  model.  Also,  it  has  been  found  that  transient  loads  may  be  measured  on  a 
dynamic  model  and  correlated  with  analytical  results  with  good  agreement. 
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Figure  2.  Complete  Model  and  Cutaway  View 
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Figure  3-  Planform  of  Original  Stabilizer  Design 
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Figure  4.  Planform  of  Final  Stabilizer  Design 
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O  MEASURED  STIFFNESS  POINTS 


Figure  5 •  Comparison  of  Fuselage  Model  Stiffness  Distribution  with  Target 

Stiffness  Distribution 
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Figure  6.  Angular  Influence  Coefficients  for  Stabilizer  Compared  with 
Calculated  Values  in  the  Chordwise  Direction 
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Figure  7.  Angular  Influence  Coefficients  for  Stabilizer  Caspared  with 
Calculated  Values  in  the  S panwise  Direction 
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Figure  8. 


Comparison  of  Calculated  and  Measured  Flexibility  Matrices  for  the 
Stabilizer  Model  Root  Springs 
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Figure  9*  Stabilizer  Node  Lines  and  Frequencies  for  Cantilever  Hoot  Restraints 
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Figure  10. 


CVC  High  Speed  Wind  Tunnel,  Test  Section  Density  Variation  as  a 
Function  of  Mach  Number 
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Figure  11.  Typical  Data  Trace  for  Elevator  Inputs  (M  *  3.0) 


Figure  13.  Time  Sequence  for  Store  Trajectory  at  M  -  3.0 
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Figure  lh.  Calculated  Rigid  Body  Mode  Shapes  for  the  Fuselage  Model  as 

Installed  on  Tunnel  Sting 
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MODEL  FUSELAGE  STA.  (IN.  FROM  NOSE) 


Figure  14. A  Calculated  Elastic  Mode  Shapes  for  the  Fuselage  Model  as 
Installed  on  Tunnel  Sting  (Continued) 
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(This  page  is  UNCLASSIFIED) 


Figure  15.  Stabilizer  Node  Lines  and  Frequencies  for  Design  Hoot  Restraints 
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Figure  16.  Comparison  of  Calculated  and  Measured  Bending  Moments  for  a  Zero 

Airspeed  Store  Ejection 
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Figure  17.  Comparison  of  Calculated,  and  Measured  Bending  Moments  for  Store 

Ejection  Input 
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Figure  18.  Comparison  of  Calculated  and  Measured  Bending  Momenta  for  Elevator  Input 


CONFIDENTIAL 


Table  1.  Ccmpariaion  of  the  Calculated  "Free-Free"  and  the  2al< 
Inatalled-On-The-Sting  Frequenclea  of  the  Complete  I’^eaU, 

Stabilizer  Model 


FREE-FREE 

FREQUENCIES 

CPS 

STING  MOUnt  1  >:.C 

FREQUENCIES  (CAL 
CPS 

1 

0 

12  67 

2 

0 

22  57 

3 

100.46 

101.03 

4 

222.57 

222.76 

5 

313.17 

3 13.68 

Table  2. 

- -  wn 

Summary  of  Complete  Model  Inertial  Data 

MODEL 

DESIGN 


MASS 

(GM) 

C.G. 

(IN  FROM  NOSE) 

'CG 

GM  IN' 

3250 

30.3 

3.25  X  10 

3124 

29.47 

3.22  X  10 

MASS  FREQ  i  ST'."  FNi 

RATIO  SATIO  i  RA  T  O 


RATIO  .SATIO 


1.03  !  0.606  I 


1.00  j^0.548 


Table  3.  HSWT  Test  Points 


RUN  NO. 

'  1 

MACH  NO.  | 

1 

3.0 

2 

3.0 

3 

3.0 

4 

3.2 

5 

2.8 

6 

2.8 

7 

2.8 

8 

3.2 

9 

2.6 

EJECTION  DATA 
Q  (PSI) 
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